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Abstract. Semiconductor microcavities, in which strong coupling of excitons 
to confined photon modes leads to the formation of exciton-polariton modes, 
have increasingly become a focus for the study of spontaneous coherence, lasing, 
and condensation in solid state systems. This review discusses the significant 
experimental progress to date, the phenomena associated with coherence which 
have been observed, and also discusses in some detail the different theoretical 
models that have been used to study such systems. We consider both the case 
of non-resonant pumping, in which coherence may spontaneously arise, and the 
related topics of resonant pumping, and the optical parametric oscillator. 



1. Introduction 

Semiconductor microcavities have been designed to greatly enhance the matter-light 
interaction strength by confining light. The confined light couples to excitonic 
resonances in the medium inside the microcavity; when this exciton-photon coupling 
exceeds the exciton and photon damping rates one finds spectrally separated normal 
modes, microcavity exciton-polaritons [|[ |^. Unlike many other examples of strong 
coupling, studied in quantum optics, this review will focus on planar semiconductor 
microcavities that confine photons in only one direction, thus leading to a continuum 
of strongly coupled modes. Due to their dual matter-light nature, exciton polaritons 
can be manipulated and studied through their light component, and have an effective 
interaction through their matter component and the nonlinearity of light-matter 
coupling. Due to the continuum of modes, the behaviour of exciton polaritons may 
be related to the statistical mechanics of interacting bosons. Thus, semiconductor 
microcavities provide an ideal system in which to study the interface between quantum 
optics, strong coupling, spontaneous coherence and quantum condensation. 

A closely related area of research, although one we will not address in this 
review, is strong coupling to single excitonic resonances i.e. excitons in quantum 
dots in semiconductor microcavities. As well as experiments on quantum dots in 
planar microcavities , experiments with confinement of photons in either three or 
two spatial directions are also performed. In the first case, OD microcavities have 
been constructed from photonic crystals jj, ^ (in which in-plane confinement results 
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from localisation on a defect in the photonic crystal, and vertical confinement from 
total internal reflection), as well as micropillars |^ (where Bragg mirrors provide 
vertical confinement, and total internal reflection provides in-plane confinement). 
In the second case, wire structures were fabricated by chemical etching Lying 
between confinement in one and three spatial directions are experiments in patterned 
microcavities, where schematically a variation in the width of the microcavity provides 
a shallow in-plane trap for photon modes ^, |T^. This results in coexistence of 
OD and 2D polariton states, separated in energies: A clear polariton spectrum has 
been seen along with the quantisation induced by a box-like confinement. Polariton 
localisation due to the intrinsic photonic disorder has also been observed . 

Our review will concentrate on macroscopic collective phenomena arising from the 
interaction between these special bosonic particles. While some of these issues have 
been addressed in other contexts — such as quantum condensation in dilute atomic 
gases 13 , and coherent quantum optics of lasers ||l^ — the combination of effects 
seen in microcavity polaritons calls for new approaches. In part, the new theoretical 
challenges arise from description of features of semiconductor microcavities, such as 
the disorder and decoherence that arise in solids, the spin structure of polaritons, 
and the effect of pumping and decay. In addition, polariton systems provide new 
opportunities for experimental probes and observations that differ from those possible 
in other systems. 

The recent experimental and theoretical research in this field can be divided into 
two main directions: Firstly, experiments which use resonant (coherent) pumping have 
been motivated by the search for all-optical ultrafast switches and amplifiers. The 
other direction is that of non-resonantly pumped microcavities, where experiments 
have pursued the search for Bose-Einstein condensation (BEG), polariton lasing and 
macroscopic phase coherence phenomena. 

A number of reviews have been already written on the subject of microcavity 
polaritons: In Refs. Jl5ljlq|, linear properties of microcavities have been analysed in 
great detail. Refs. ]17|, |18{ review problems of non-linear optics and the theoretical 
framework necessary to study resonant pumping, parametric amplification and 
oscillation. In addition, two books ^ and two special issues ^ have also 
been published. In this review, we will focus attention on those issues of modelling 
microcavity polaritons that are of particular importance in understanding spontaneous 
coherence and condensation in such systems. We will address the relation between the 
various theoretical approaches that have been used, and discuss the limits in which 
they become equivalent. We will also discuss in some detail the relation between 
coherence, condensation, lasing and superfluidity in experimental systems which are 
finite, two-dimensional, decaying, and interacting, and thus differ from the Bose- 
Einstein condensation of ideal three-dimensional bosons. Finally, we will discuss the 
relation of these features both to the resonantly pumped polariton system, and also 
to other experimental systems in which similar issues of coherence and condensation 
in complex systems are addressed. 

The review is divided overall into non-resonant pumping in Sec. |^ and resonant 
pumping in Sec. |^. Within Sec. |^, we first present the experimental development 



of the subject in Sec. 2.1, and then discuss the theoretical approach, dividing our 



discussion into the question of choice of model in Sec. p.2|, choice of treatment (i.e 



thermal equilibrium, rate equations, etc.) in Sec. 2.3, and the phenomena predicted 
(i.e. experimental signatures, conditions for condensation) in Sec. 2.4. Within Sec. || 
we again divide into a summary of experimental progress in Sec. 3.1, and a discussion 




Figure 1. Schematic diagram of a microcavity, formed by a pair of distributed 
Bragg reflector stacks, with quantum wells at the antinodes of the cavity photon 
mode. 



of the additional theoretical issues relevant only to the resonantly pumped case in 
Sec. 3.2, Section|^ finally draws comparisons to phenomena seen in other experimental 
systems, and briefly summarises our discussion. 



1.1. Introduction to microcavity polaritons 

Before discussing experiments and theories of coherence in microcavity polaritons, 
we provide here a brief introduction to the systems considered, and to microcavity 
polariton modes. Fuller introductions can be found elsewhere |l^, The 
semiconductor microcavities we discuss are constructed from distributed Bragg 
reflectors, containing alternating quarter wavelength thick layers of dielectrics with 
differing refractive indices. Due to these Bragg reflectors, the cavity contains a 
standing wave pattern of confined radiation. As illustrated in Fig. |l|, quantum wells 
(QWs) are placed at the antinodes of this standing wave, thus maximising the coupling 
between photons and excitons confined to the quantum wells. 

Because the photon modes are confined to the cavity, the volume associated with 
the radiation mode is small, and so the exciton-photon coupling is strong. This strong 
coupling means that rather than considering the exciton-photon coupling as leading 
to radiative decay of the excitons, the exciton and photon modes are instead mixed, 
to form new normal modes: lower and upper polaritons. At the simplest level, one 
can write the exciton-photon Hamiltonian in terms of operators V'k creating photons 
and Z?Jj. creating excitons, with k labelling the 2D in-plane momentum. Thus: 

r )(X)' 

[We have set h = 1 here and throughout.] Here, Wk is the energy of the photon mode 
confined in the cavity of width w, giving: Wk = (c/n)-\/fc^ -f- {2ttN/w)^, with n is 
the refractive index, and N the index of the transverse mode in the cavity. For the 
situation in Fig. |], = 2. For small fc, the energy can be written as Wk = tUQ + k'^ /2m, 
where m is an effective photon mass m — {n / c){2'k / w) . In the absence of disorder. 
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Figure 2. Schematic polariton spectrum. Left as a function of emission angle, 
9 = sin~^(cfc/a;o); Right, as a function of momentum on a logarithmic scale, 
showing the full exciton dispersion. (Plotted for M = O.OSrrie, m = 3 X 10~^me, 
= 26meV, 5 = 5.4meV and ujq = 1.7eV.) 



the exciton energy in the QW is ek = £o + fc^/2Af, where M is the total exciton 
mass, and eo = Ecv — Ti-ex comes from the conduction- valence band gap Ecv including 
QW confinement and the exciton binding energy (Rydberg) Ti-ex- For convenience, 
we define the bottom of the exciton band, Eq a-s the zero of energies; and denote the 
detuning between exciton and photon bands as (5 = t^o ~ ^o- Finally, the off diagonal 
term describes the exciton-photon coupling, where is the Rabi frequency. 

Then, diagonalising the quadratic form in Eq. (|^) gives the polariton spectrum: 



=2 




(2) 



This spectrum is illustrated in Fig. ||. It is shown there both as a function of 
momentum fc, and also as a function of angle. The angle corresponds to the angle 
of emission of a photon out of the cavity; since the in-plane momentum and photon 
frequency are both conserved as photons escape through the Bragg mirrors, one may 
write (eo -|- E^^) sin[9) — ck, which (since typical values of k satisfy wq, £o ^ ck) can 
be approximated as loq sin(0) = fc^/2m. 



2. Non-resonant pumping 
2.1. Summary of experiments 

Optical properties of semiconductor microcavities have been the subject of extensive 
experimental investigations since the first observation of the strong coupling regime 
by Weisbuch et al. Much of the experimental research has concentrated on 

III-V materials, mainly GaAs/AlGaAs structures, or on II- VI materials, such as 
CdTe/CdMnTe/CdMgTe structures. The main aim of the experiments described 
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here as non-resonant pumping has been to start with incoherently injected polaritons, 
and observe spontaneous coherent processes emerging from incoherent injection of 
polaritons: polariton degeneracy, final state stimulation, and ultimately polariton 
BEG. 

The authors of the earliest report of non-linear emission in the presence of 
strong coupling in GaAs microcavities, suggesting final state stimulation characteristic 
for bosonic particles, later withdrew those conclusions | p5[ as further experiments 
showed that the threshold for non-linear emission occurred in that case after the 
crossover to the weak-coupling regime; thus the non-linear emission should have been 
attributed to photon lasing. The first unambiguous observation of polariton bosonic 
stimulation was in CdTe microcavities ||2^ consisting of 16 quantum wells with a Rabi 
splitting of around 23meV. Two distinct stimulation thresholds were observed with 
increasing intensity of continuous wave pumping, as shown in Fig. ^. As the pumping 
intensity was increased above the first threshold, non-linear emission at energies close 
to the bottom of the lower polariton branch was clearly seen. The second threshold, 
reported for much higher intensities, was connected with a weak coupling electron- 
hole lasing mechanism. Further investigation |2^, |2|] showed that the first nonlinear 
threshold — in the strong coupling regime — was due to stimulated scattering to 
the ground state. Very shortly after publication of Ref. , non-linear emission was 
seen in a single QW GaAs microcavity characterised by 3.5meV Rabi splitting. 
However further investigation jSOt showed that this nonlinearity had emission varying 
as the square of pumping intensity, and a threshold that occurred for occupation 
factors much less than one, and so this nonlinearity was associated with increase of 
exciton-exciton scattering, rather than final state stimulation. Final state stimulation 
in III-V materials was demonstrated by pump-probe experiments in Refs. |3l[ |. 

The stimulated scattering to the ground state, and non-linear build-up of 
lower polariton population was the first step towards demonstration of spontaneous 
coherence and thermalisation — characteristic of quantum condensation. However, 
the big challenge to realising a condensed polariton phase was the finite (though very 
large) quality of the cavity mirrors, and the resultant short polariton lifetime, of the 
order of picoseconds. In addition, due to the 'bottleneck effect', ||3|] the relaxation 
of polaritons to the zero momentum state was delayed, hindering the creation of a 
thermal population in the lowest energy states. The first investigation of the coherence 
properties of emitted light above the threshold for non-linear emission in strong 
coupling was based on the measurement |]33| , Q of the second order coherence function, 
g2{t = 0), which would take a value of 52(0) = 2 for a thermal state, and 52(0) = 1 
for a coherent state js^. A decrease of 52(0) from 1.8 to 1.4 as pumping power was 
increased from threshold to 20 times threshold power was seen in a system of 12 GaAs 
quantum wells placed at the antinodes of light in a GaAs/AlGaAs microcavity, giving 
14.9meV Rabi splitting. This was followed by a report of a characteristic change in the 
momentum space distribution above threshold ||3^, as shown in Fig. ^, and a blueshift 
of the polariton dispersion p7[ . Time-resolved photoluminescence measurements were 
also performed for a GdTe microcavity with lO.SmeV Rabi splitting |Q under non- 
resonant pulsed excitation, which were able to monitor the buildup of a large polariton 
population in the k = state. Analysis of the time dependence showed that, below 
threshold, the dynamics of the k = polaritons follows closely the population of 
cold reservoir excitons; the relaxation from high energy exciton states resonant with 
the pump to these reservoir exciton states had a characteristic relaxation time of 
30ps. [Note that this time is significantly shorter than the 150ps observed in similar 
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Figure 3. Panel (a): Reflectivity spectrum, showing location of lower and 
upper polariton modes in the absence of pumping. Panel (b): Photoluminescence 
as pumping power is increased, a threshold for nonlinear emission is seen at 
40kW/cm-^, while the system clearly remains in the strong coupling regime. In 
panel (c), at much higher pumping powers (10,000kW/cm2), a second threshold, 
to electron-hole lasing is seen. [From Copyright (1998) by the American 
Physical Society] 

experiments [|9| with GaAs based microcavities.] Above the non-linear threshold, 
the population of the reservoir excitons was found to be clamped, and the polariton 
relaxation dynamics became faster, with the maximum of polariton emission at 70ps 
delay after the initial pulse. This delay was further decreased for higher excitations 
powers. Together these provide evidence of stimulated exciton-exciton scattering to 
the lower polariton states. 

The first evidence of spontaneous first-order coherence in an incoherently pumped 
microcavity was seen in a 16 QW CdTe microcavity with 26meV Rabi splitting |Q 
under non-resonant pulsed pumping. An interesting feature of this particular 
experiment was that the non-linear emission was at k 7^ and so resulted in an 
emission ring at an angle of around 17°; this was associated with the small size of the 
excitation spot (3 /im). (Note that in later experiments with larger excitation spots 
on the same sample condensation was at k = [0.) The first-order coherence was 
investigated by spectroscopic imaging of the far-field emission. Two momentum space 
images were superimposed giving fringes (as a function of momentum k) with over 75% 
contrast above threshold and up to 35% below threshold. In a later publication of the 
same group experiments on a 4 QW CdTe microcavity characterised by 13.2meV 
Rabi splitting showed macroscopic occupation of the k = state characterised by 
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Figure 4. Momentum space distribution of lower polaritons above (main figure) 
and below (inset) the nonlinear threshold. [From Ref. Copyright 2003 

National Academy of Sciences] 



narrowing above threshold of the polariton emission Hne to a Hnewidth below that 
of the cavity photon mode. Near-field images showed modulation of the polariton 
spatial distribution, revealing the effect of photonic disorder. The next challenge in 
the search for spontaneous condensation was to see similar effects, but accompanied 
by a thermal distribution of k polaritons. 

Due to the short polariton lifetime and the 'bottleneck effect', |Q the realisation 
of equilibrium population has proven to be challenging. Progress came from 
observing ^ that thermalisation processes due to particle-particle scattering 
can be dramatically increased both by increasing the value of the (non-resonant) 
pump power, and also by positively detuning the cavity energy above the excitonic 
energy. Large positive detuning makes polaritons more excitonic and increases their 
scattering rate. Time- and angle-resolved spectroscopy on a sample consisting of 12 
GaAs QWs characterised by 14.4meV Rabi splitting |Q showed that for positively 
detuned cases, where the thermalisation time increases while decay time decreases, 
the thermalisation time can reach around one tenth of the polariton lifetime and that 
lower polaritons remain in thermal equilibrium with the phonon bath for a period of 
about 20ps. Finally, a comprehensive set of experiments showing clear evidence for 
condensation of cavity polaritons was performed in a CdTe |^ structure consisting of 
16 quantum wells giving 26meV Rabi splitting. Above the threshold pumping density 
they observed: a massive occupation of the k = mode developing from a polariton 
gas in thermal equilibrium at 19K (shown in Fig. an increase of temporal coherence 
from 1.5ps below threshold to up to 6ps above threshold; the build-up of long-range 
spatial coherence over the whole system size with contrast of interference fringes 
from less than 5% below threshold to 45% above threshold; linear polarisation of the 
emission. Linear polarisation had been predicted to appear in the condensed state |Q| , 
and its appearance gives evidence for the single state nature of the condensate, though 
since the direction of this polarisation was pinned to a crystallographic direction the 
polarisation direction symmetry was not spontaneously broken. Evidence that the 
polarisation of light is pinned to one of the crystallographic axes, independently of 
the excitation polarisation, was also independently observed in experiments on both 
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CdTe and InGaAs microcavities ji^l ■ These results were ascribed to birefringence 
in the mirrors and cavity. Shortly following the work in Re f. |41[ , a similar non- linear 
build up accompanied by a linear polarisation was seen |[47|] in GaAs structures in 
stress- induced traps Q. Also the second order coherence function, 52 = 0), has 
been measured in the CdTe structures studied in Ref. |41 where, in contrast to 
observations reported in Ref. |3^, 32(0) was found to be around 1 at threshold and 
then increased up to around 1.4 at powers 10 times threshold. This effect has been 
attributed to the phase diffusion due to interactions. 
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Figure 5. Polariton occupation vs energy, sliowing evolution from below to above 
threshold power, with the distribution remaining thermalised at a temperature of 
19K [From Ref. 0] 



Wide-band-gap semiconductor structures based on group-Ill nitrides, such as 
GaN based cavities have recently attracted considerable interest (see, e.g., |5^, 
|3|,|54[||, Hl^). The main advantage of these structures over II- VI and other III-V 
materials lies in the large exciton binding energy (around 26meV for bulk structures, 
and over 40meV for narrow quantum wells) and the large coupling to the photon field, 
which makes them ideal systems for the realisation of functional devices operating 
at room temperature. Although the study of cavity polaritons in group-Ill nitrides 
microcavities is still in its infancy, strong exciton-photon coupling in a bulk GaN 
cavity js^, ^ |5l[ and in a quantum well cavity have been reported. In both cases, 
the substantial inhomogeneous broadening of the excitonic and photonic lines play a 
key role in establishing the conditions for reaching strong coupling. Very recently, a 
non-linear build-up of polariton emission accompanied by an increase in the first-order 
temporal coherence, and a spontaneously chosen linear polarisation (independent of 
the apparatus and different between measurements) has been reported to occur at 
room temperature in GaN bulk microcavities |p3| . 
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2.2. Theoretical models 

In this section, we will discuss the different models that have been used to describe 
microcavity polaritons and study their condensation. We wish to separate clearly two 
aspects of theoretical description of polaritons; the first aspect is the choice of model, 
the subject of this section, the second aspect is how that model is treated, which 



will instead be covered in Sec. 2.2. After having addressed these points, we then in 



Sec. 2.4 discuss the various theoretical predictions of conditions for condensation, and 



of possible signatures. Readers who are not interested in the details of how the system 



is theoretically modelled should jump to Sec. 2.4. A model starting from electrons 
and holes, taking into account their Coulomb interaction to form bound excitons, 
their coupling to light, and the effects of disorder would describe polariton systems 
exactly, but is too complicated to allow any clear understanding of the important 
features associated with condensation to be gained. Therefore, it is appropriate to use 
simplified models, that exaggerate some features of the real system, and neglect others. 
In judging which model is appropriate to address a given problem, it is important to 
understand how the model relates to the underlying microscopic model of electrons 
and holes, and so we shall start by discussing this microscopic model. 

2.2.1. Microscopic electron-hole Hamiltonian In this section, we discuss the 
underlying description of microcavity polaritons formed from photons confined to 
a two-dimensional cavity, interacting with electrons and holes in two-dimensional 
quantum wells ^ . 

H Hq\y -f^coul ^" -^disorder + ^fphoton + ^^di pole (3) 

Consider first the electrons and holes, we have 

= 5]ekCkCk + ek4«k (4) 



iJdisordor- ( dv[We{r)c\v)c{v)-Wh{v)v{r)v^{v)\ . (6) 



coul 



Here cj^ (wj^.) create electrons in the conduction (valence) bands, which have 
dispersions, (Sk)- Since the "empty" state is a filled valence band, it is more 
convenient to describe the valence band via the operator Wk which creates a hole — 
i.e. a missing electron. The density of electrons (holes) is given by Pq — ^j^. c\^^^c^ 

(Pq — X^k ''^k^k+q)- "^^^ factor I /A, where A is the quantisation area of the cavity, 
appears explicitly because the Hamiltonian has been written as a sum over momentum 
labels; this factor plays no role in any final answer, and is absorbed in the definition of 
dk if summation is replaced by integration. Note also that in general there should be 
a dependence on the electron and hole spin degrees of freedom, that we neglect here. 
The last term, Eq. (||), describes the disorder potential acting on electrons and holes, 
e.g. due to well-width fiuctuations and alloy disorder. In general, disorder can act 
differently on electrons and holes; in practice for the materials used, the energy scale 
of disorder is less than the binding energy, so disorder does not dissociate excitons |^ . 
If the exciton binding energy is significantly larger than a characteristic energy scale 
of disorder, then as described in Ref. BCfl one can factorise the wavefunction into a 
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centre of mass wavefunction, and a wavefunction of relative electron-hole separation. 
Then, in the equation for the centre of mass wavefunction, one has an effective disorder 
potential that is the result of convolving the original disorder with the wavefunction for 
relative electron-hole separation. As a result of this convolution, the effective disorder 
potential as seen by the exciton centre of mass wavefunction is smoothed over the 
scale of the exciton Bohr radius pH |. 

Turning now to the interaction with the photons, 

i/photon = E^q^^q + / *W^ph(r)V^(r)V'(r) (7) 
q •' 

q,k ' 

In Eq. , the quantisation volume for the electromagnetic field has been factored into 
ALw, where Lw is the width of the cavity, and A the quantisation area as discussed 
above. The term /^cu is the inter-band dipole matrix element, which can be calculated 
given the Bloch wavefunctions of the two bands. The term Wph(r) in Eq. (Q) describes 
photonic disorder, which can arise due to roughness of the Bragg mirrors — i.e. due to 
layer width fluctuations (monolayer mismatch), or crystal dislocations |ll], The 
effects of this photonic disorder, and of the exciton disorder introduced above, can 
be quite different. The photonic disorder is generally on large length scales (typically 
of the order of a micrometre), comparable to the size of the excitation spot, and 
so it is primarily associated with the spatial inhomogeneity of polaritons seen in 
experiment |42 . In contrast, as discussed in Sec. 2.2.3| , excitonic disorder is on much 



shorter length scales (typically of the order often nanometres for CdTe), and thus does 
not affect the spatial polariton density profile; however excitonic disorder does have a 
significant impact on the distribution of excitonic oscillator strengths. Although the 
excitons are localised, in the absence of photonic disorder the polaritons formed consist 
of a superposition of many different localised excitons and extended photon states, 
and thus one may form delocalised polaritons from localised excitons [Q. We will 
not explicitly discuss the effects of pho tonic disorder further, however the discussion 



of condensation in a trap in Sec. 2.3.3 can apply also to trapping in disorder, as well 
as any deliberately engineered trapping. 

The Hamiltonian in Eq. (y|) already contains a number of important 
approximations, which should be discussed. The interaction of photons with electrons 
and holes makes use of both the dipole approximation, and the rotating wave 
approximation [^3| , Chapter 10]. The interaction strength here is written in the dipole 
(length) gauge. The choice between the dipole (length) gauge and the Coulomb 
(velocity) gauge is not arbitrary, as the terms assigned as describing free particles 
(without interaction with radiation) are different in each gauge |6^. This 

point is worth stressing, as the electromagnetic interaction between excitons is split 
between the direct Coulomb term, and a photon mediated term. Thus the choice 
of gauge affects also the Coulomb interaction [Eq. (||)] , controlling which parts of it 
are absorbed into the definition of exciton states, which parts are associated with the 
"photon" operators — in the Dipole gauge, the fields V'q are quantised modes of the 
electric displacement — and which should be written as some effective exciton-exciton 
interaction ]6^ . The relation between Coulomb interaction and photon mediated 
interaction is complicated here because the resonant photons are confined by the 
DBR (distributed Bragg reflector) mirrors, while the static Coulomb term is modified 
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much less strongly by the mirrors. When one comes to exciton states, it is therefore 
important to be aware that the choice of gauge affects both the exciton-photon 
coupling strength, and the form of the inter-exciton Coulomb interaction, and that 
these two are not separate. 

In the next two sections, we will discuss the main two classes of effective 
Hamiltonians, derived from this full Hamiltonian, used to study microcavity 
polaritons. The differences between these effective Hamiltonians can be seen as 
the result of regarding different terms as important; i.e. which terms are treated 
exactly, and which perturbatively. In both cases, the first step involves changing 
from electrons and holes to bound excitons — i.e. solving the wavefunction for the 
relative coordinates. The differences then arise from considering in one case next the 
effect of disorder, giving localised states, and then approximating the inter-exciton 
Coulomb term by exclusion — this leads to the boson- fermion model discussed in 



Sec. 2.2.3 — or alternatively, treating the Coulomb term via a quartic exciton-exciton 



interaction term, then coupling to light, and then treating disorder perturbatively or 



not at all — this leads to the weakly interacting boson model, discussed in Sec. 2.2.2 . 

As will be discussed further below, in the low density limit, many features of 
these models are similar. However, the different models emphasise different features: 
The boson model can effectively describe the case where the dominant interactions are 
exciton-exciton Coulomb interactions, while the boson-fermion model instead has the 
saturation of the exciton-photon coupling as the dominant interaction. As such, these 
different models may be appropriate in different contexts. For example, to describe 
the lower polariton blue-shift, and comparable upper polariton red-shift seen, e.g. 



in Refs. |41, the effects of the saturation interaction are required. Further, the 
different models have been developed in different directions, for example the effects of 
exciton spin, and thus the polarisation dynamics have so far only been considered in 
the bosonic model (see e.g. Refs. ||6^, Q and Refs. therein). 

2.2.2. Weakly interacting boson models A weakly interacting Bose gas model of 
polaritons can be achieved by making an Usui transformation p8| , |69| , choosing 
the bosonic operators to represent bound exciton states, and then truncating the 
interaction terms at fourth order ||7^. This results in an effective Hamiltonian 
describing bosonic excitons coupled to photon modes: 

n, 



k 

n 



- E K'-q^+q^kV'k'+^L^qi^Uq^k^k' 



k.k'.q 



+ E ^^<c^L-.^k'^k, (9) 



k,k' 



Here D^^ creates a bound exciton of energy £k, "^k creates a cavity photon of energy 
LOk, and flfi is the effective exciton-photon coupling strength, or Rabi splitting. By 
measuring energies from the bottom of the exciton dispersion, we may write — 
k'^ /2M, and expanding the photon dispersion to quadratic order in momentum, one 
may approximate — k'^ /2m + S, with S the exciton-photon detuning. The quartic 
terms in Eq. (^) are divided into exciton-exciton interactions, ?7k-k',q, the strength 
of which can also be found by calculation of the Coulomb exchange term |7^, in 
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the Born approximation, and a "saturation term" (second line), which decreases the 
exciton-photon coupHng at large exciton densities due to the fermionic character of the 
excitons [ [70| . These quartic terms arising from the Usui transformation can be seen 
as an expansion of the underlying fermionic operators in powers of bosonic operators; 
this expansion is controlled by the small parameter of the number of excitons per Bohr 
radius. Note that in general these terms depends also on the spin degrees of freedom 
of the constituent electron and holes. For a derivation of the dependence on spin of 
the Coulomb terms, see, e.g., Refs. |73, [74|. 



This approach takes into account the intra-exciton Coulomb term, in forming 
bound excitons, and the inter-exciton Coulomb terms as an effective quartic 
interaction. The Hamiltonian in Eq. however neglects disorder acting on the 
exciton states, and as a result finds that each exciton state couples to a sin gle ph oton 



state, with conserved momentum. However, as discussed below in Sec. 2.2.3 , and 
in Refs. ||7^, exciton disorder will modify this picture. Including disorder, one 
finds a distribution of energies, and at each energy a distribution of exciton-photon 
coupling strengths. The exciton states which have the largest coupling strength to the 
low momentum photons are found to be at energies just below the exciton dispersion 
edge. Although these state s are not the most localised, i.e. are not the states far in the 
Lifshitz tail (see later Sec. 2.2.3 ), they are below the band edge and therefore are still 
quite strongly localised, they decay quickly at long distances, and they result from 
exciton wavefunctions concentrated around minima of the potential. The low energy 
polariton modes will be formed from a superposition of many such localised exciton 
states. 

In addition, the saturation term in this model, which describes the reduction of 
exciton-photon coupling is taken only to the lowest order. This is sufficient at low 
enough densities; however as discussed more fully in Sec. 2.2.3, including effects of 
disorder the density at which these saturation effects become important can be much 
lower than the Mott density. Thus, a quartic description of saturation may become 
inadequate at modest densities, close to those already studied experimentally. In 
addition, most bosonic models of polaritons further simplify Eq. (^, replacing the 
momentum dependent interaction C/k-k',q with its strength at k = k',q = 0. This 
strength is the interaction between two excitons in the same single particle momentum 
eigenstate. If exciton eigenstates are localised, it is not obvious that replacing 
all exciton-exciton Coulomb interactions with an average strength (calculated from 
delocalised exciton wavefunctions) is appropriate. In addition, the dominant Coulomb 
interaction between localised, and therefore non-overlapping, exciton states may well 
be due to the direct dipole-dipole interaction, rather th an exc hange terms (as it is in 
the clean case). The boson- fermion model discussed Sec. 2.2.3 handles this interaction 
differently — it includes strong on-site repulsion, and neglects inter-site repulsion; this 
limit is clearly also an exaggeration, and the true effects of Coulomb will be between 
these two extremes. 

It is worth noting parenthetically that a constraint on exciton density 
X]k(-^k-^k) ^ Pso-t required to make the Hamiltonian in Eq. (||) stable. Without 
such a constraint the free energy is unbounded from below, i.e. for \^) = 

exp ^AtAq + PO}^ |0), the free energy F = {H — ^N) corresponding to Eq. (^ is 



F = {d- ^i)\x\' - + nn^ixn i ~ ^ + (lo) 

V Psat / 2 

The minimum free energy can be found for real A, (3, and so re-parameterising these 
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as A = a;sin(x), j3 ■ 



xcos{x) 
Uo 



4 ? 

X COS X 



the quartic term in Eq. (|T^) goes like: 



1 + cos(2x) 



■ sin(2x) 



(11) 



UoPsat 

For any non- vanishing f2fl, there is a value of x lor which this is negative and so 
unstable. Physically this instability is cured by restoring higher order contributions 
of the saturation interaction which prevent X]k(-^k^k) P^at- Practically the above 
instability can be avoided if one diagonalises the quadratic part of Eq. (j^) , and then 
projects onto the basis of lower polariton states iTn]. By writing: 



cos 
sin ^k 



- sm Wk 
cos 0k 



4 



(12) 



here LI , Ul create lower and upper polaritons respectively, and cos^k, sin^k are the 



standard Hopfield coefficients 
Eq. (^, one must choose 

tan(26'k 



t^k - £k 

with o^k and £k as defined following Eq. 



|59f|. In order to diagonalise the quadratic part of 



(13) 



Having diagonaliscd the quadratic part, 



one may project onto the lower polariton basis for the quartic part, giving the effective 
lower polariton Hamiltonian: 



LP 



Eh" 



LP rt r 
^k^k 



E 

k,k',q 



off 

k,k',q^k 



L 



q-^ic'-q-^k'^k 



(14) 



(wk + Sk) - V {uJk - ek 



(15) 



off 

k,k',q 



■ COS 6'k+q COS 6'k [cos ^k'-q sin 6'k' + sin 6'k'-q cos 6'k' 



'2psat 

— COS 9k+q COS ^k COS 0k'-q COS 6^k' 



(16) 



Note that in order for the neglect of upper polaritons to be valid, one must be 
at temperatures significantly smaller than the Rabi splitting. This requirement of 
temperature can be translated to a requirement of low densities if one is interested 
in phase transitions: the density must be low enough that the Bose condensation 
temperature at that density is much less than the Rabi splitting. It can be shown |Q 
that this latter requirement means one should have fewer than one polariton per 
wavelength of light; such a density is already exceeded in current experiments. 

The Hamiltonian ( |l^ has an effective k dependent interaction strength due to 
the change of Hopfield coefficient along the lower polariton branch — i.e. Coulomb 
interaction becomes stronger as the polariton becomes more excitonic, and saturation 
interaction is strongest nearest to equal photon and exciton components. Preserving 
a k dependent coupling strength requires one to think carefully about regularisation. 
In atomic gases, the weakly interacting Bose gas model is generally studied with a 
contact interaction, as is appropriate when the scattering length is much less than 
the de Broglie wavelength; this is renormalised by matching the scattering length to 
the experimentally measured quantity [ p^ . If the interaction is instead found from a 
microscopic theory, as it is the case here, that microscopic theory must also describe the 
regularisation of the interaction at large momentum, as a single measured scattering 
length would not allow fitting of the different momentum dependencies associated 
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with Coulomb and saturation terms. In practice, this means any attempt to preserve 
the effect of Hopfield coefficients on the interaction must also take into account the 
decrease of both Coulomb and saturation effects for large exchanged momenta. 

The limits of validity of this Hamiltonian come from several sources; the 
requirement for density to be less than the Mott density l/a\ is the easiest to 
understand, but also the most easily satisfied. Neglect of the upper polariton required 
temperatures less than the Rabi splitting (which translates to densities less than one 
polariton per square wavelength of light) ; but inclusion of the upper polariton leads 
to instabilities, which would require higher order terms in the Hamiltonian to restore 
stability. Thus, consideration of the phase boundary at high densities, at which 
the naive estimate of the transition temperature would be comparable to the Rabi 
splitting, would require a treatment beyond that considered in this section. Thus, in 
the next section we discuss an alternative model that should be valid at these higher, 
yet still experimentally accessible densities, and also takes account of the effect of 
disorder on the saturation interaction. 



2.2.3. Boson-fermion, and generalised Dicke model By considering first the effects 
of disorder acting on the excitons, one finds that in 2D systems the effect of disorder is 
particularly profound and that formally any arbitrarily small amount of disorder leads 
to localisation |7^, However, the character of the states changes significantly with 
energy. At high energies states may be described as a random superposition of plane 
waves with the same modulus of momentum, and localisation effects arc weak. At 
very low energies, well below the band edge, the Lifshitz tail states [^0[ ^ have 
a nodeless form, localised in deep minima. The changing nature of the exciton states 
with energy also changes their oscillator strength and the exciton states that 

couple most strongly to the long wavelength radiation modes are those just below the 
band edge, for which localisation effects are important. As a result, those exciton 
states which contribute most to the relevant (thermally populated) polariton states 
are effectively localised exciton states ^ . 

This localisation may also be expected to modify details of the inter-exciton 
Coulomb interaction term compared to the clean picture j?^, ^ . Considering strongly 
localised exciton states, since exchange requires wavefunction overlap, one expects 
a difference between the strength of on-site Coulomb repulsion — i.e. interaction 
of excitons localised in the same potential fluctuation — as compared to inter-site 
interactions. Taking the extreme form of this difference — i.e. on-site exclusion and 
neglect (or perturbative treatment) of the inter-site interaction — leads one to a 



generalisation [[75|, [76| of the Dicke model |85[ g6|, |87|, describing two-level systems 
coupled to a bosonic field: 

H=J2 e-^S + E ^P^P^P + E E i9'^-^P%S+ + H.c.) . (17) 

a p ^ a p 

Here Sq, is a spin 1/2, representing a two-level system, where ||q) is the ground 
state — i.e. no exciton on site a — and jta) indicates the presence of an exciton 
on site a. Such a model has also been studied in the related context of spontaneous 
superradiance |9[ |9^ ; it was however later shown that including higher order 
terms beyond the dipole approximation prevent the superradiant transition of the 
vacuum state of such a model. No such problem however occurs when one considers 
the system in contact with a reservoir that fixes particle density — the effects discussed 
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in Ref. apply to the stability of the vacuum state, i.e. with chemical potential going 
to negative infinity. 

It is often convenient to represent the two-level systems as two fermionic states 
so that the ground state is \[a) = |0)i ^^id the excitonic state Ita) = |0) = 
h\aa |g-S.). Imposing a constraint on total fermion occupancy, h]j}a + ajjCQ. — 1, 
eliminates the unphysical states |0) and a^J)]^ |0), thus giving the Hamiltonian 

^ = XI Y (^"^° + °^""") + Xl'^p^p^p 

Q p 

nY^Y. [ffc,p^p6iaa + H.c] . (18) 

^ a p 

This formalism is easy to use, as one may show that the constraint preventing 
double occupation can be easily incorporated in the imaginary time path integral 
formalism by shifting the fermionic Matsubara frequencies according to 

e„ = (2n + l)7r//3 ^ en ^ (2n + 3/2)7r//3 . (19) 

It is important not to confuse these fermionic states (which represent the two levels of a 
two-level system) with the conduction and valence band states in Eq. (^. While b]^a^ 
creates an exciton, one should not think of (a,^) as creating an electron (hole) — 
i.e. one cannot write foj, [a^) as a linear combination of the electron creation operators 
cjkJ '''k ('^k' "^k) '^^ Eq. (^. If 6jj (ttg,) were a linear combination of electron (hole) creation 
operators, then h]^a^ would create an electron-hole pair, but without any correlation 
between the position of electron and hole — i.e. without excitonic binding. Instead, 
the relation between b\^a^ — the fermionic representation of saturable excitons — and 
the underlying electrons and holes is as discussed later in Eq. ( p^ and Eq. (p9[). 

This model naturally allows one to consider a distribution of excitonic energies, 
and a distribution of excitonic oscillator strengths for each given energy, which are set 
by disorder To perform such a calculation, Eq, and ga,-p should be calculated 

by solving the Schrodinger equation for the exciton centre of mass coordinates in a 
random disorder potential, 



V 



2 



$„(R) =£„$„(R). (20) 
and then calculating the oscillator strength from 



ffa.p = eAic^-/- ^</?ls(0)$a,p. (21) 

Energies and oscillator strengths calculated this way are shown in Fig. |[ It is of 
course perfectly possible to write a theory of disorder-localised exciton states that are 
bosonic modes — and as will be shown below, a bosonic model can be extracted from 
this boson-fermion model at very low densities — however in such a treatment it is 
important to consider, as discussed above, how the change of exciton wavefunctions 
modifies the Coulomb interaction between excitons. 

A relationship between the model of this section, and that of the previous section, 
may be established in the limit of low densities, by considering a Holstein-Primakoff 
transformation of the Hamiltonian in Eq. (p7|); i.e. 



Si^DiD^-]^, S+=Di^l-DiD^, S-^{S+)^ (22) 
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Figure 6. Distribution of energies, and energy dependence of exciton-light 
coupling strength, for excitons in the presence of disorder. Black points mark 
the energy and coupling strength of individual exciton states (combined from 
160 separate disorder realisations); green lines show density of states, red lines 
show mean square coupling strength. Panel (a) [From [P5| Copyright (2006) by 
the Anwican Physical Society] coupling to zero-momentum photons. Panel (b) 
[From |76|.] coupling to high momentum photon states (much beyond relevant 
states). Inset illustrates comparison of numerical values for the mean squared 
oscillator strength to the first Born approximation. 



Then, assuming the occupation of excitons to be small (i.e. (-Dq-Dq) ^ 1), one may 
expand Eq. ( p7| ) to get: 

a p 

+ f 3",pV'p^l f 1 - IdID^) + H.c.) . (23) 



Comparing this to Eq. (^ shows that a bosonic model derived in this way has certain 
differences to the standard bosonic model; it obviously neglects the inter-exciton 
Coulomb term, as this was neglected in Eq. (|l7|), and the exciton energies are set 
by localised states in a disorder potential Eq, rather than k^/2M. Less obviously, but 
more importantly, the saturation interaction term is significantly stronger than would 
be suggested by Eq. (^; in that case, the mean- field energy shift at polariton density 
n is of the order of 

SE^[, ^ nnnal^. (24) 

In contrast, the term in Eq. ( [2^ ) is of the order of 

6E^P ^ nnnf, (25) 

where ^ (MWp)^^^^ is a characteristic length scale of the disorder potential. This 
energy shift is important as it relates the observed lower polariton blue-shift to the 
polariton density, and so is important in the interpretation of experiments. This result 
is valid at low temperatures; at higher temperatures one can show [|76| , that should 
be replaced by {MkeT)-'^/^. The appearance of this temperature-dependent 

length scale would not arise from a model that included only bosonic lower polaritons. 
Comparison of the equilibrium transition temperatures of the two models is discussed 



later, in Sec. 2.4.1 
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2.2.4- Comparison of models As is clear from the above discussion, the Bose- Fermi 
model in some sense encompasses a bosonic model. However, its derivation led 
naturally to the inclusion of saturation interaction, but as yet no generalisation 
including long-range Coulomb interaction has been studied. The question of 
comparing the models is therefore not so much whether one model is right or 
wrong, but whether interaction effects beyond a quartic boson-boson interaction are 



important, and so whether a description like that of Eq. (17) is necessary. At low 
enough densities and temperatures (i.e. temperatures a small fraction of the Rabi 
splitting) it is clear such a description is not necessary. However, the definition of 
"low enough" that is derived from studying when the (equilibrium) phase boundary 
of Eq. ([l^) is reproduced by a bosonic theory suggests that low enough means exciton 
separation of the order of the wavelength of light |9j, [t^, rather than, e.g. the 
exciton Bohr radius; and temperatures of the order of tenths of the Rabi splitting. 

As an alternative way to resolve the question of which approximate Hamiltonian, 
Eq. ( p^ or Eq. (^), is most appropriate for a given physical system, one can propose the 
following clear, but technically challenging approach. From both Hamiltonians, one 
can construct an approximate ground state, which can then be rewritten in terms of 
electrons, holes and photons. In both cases, we consider generalisations of the coherent 
state, which for a simple structureless boson field would be written exp(AL^) |0). 
This leads to two different trial wavefunctions for the electron-hole-photon system. 
While this is not a simple exercise — and would in fact require extensive numerical 
computation — it is a useful gedanken comparison to highlight the distinctions. Let 
us consider first the trial wavefunction appropriate to the Hamiltonian of Eq. (p^. 
Taking |0) as the filled valence band, we have 



l^Bose) = e^^« |0) ; 4 = cos(fo)^i + sin(eo) ^(^X^-q- (26) 

q 



At low densities this wavefunction has a simple interpretation; i^(q) is the bound 
exciton wavefunction, and the controls the exciton and photon fractions of the 
lower polariton; i.e. the term in brackets is the lower polariton creation operator, and 
this is a coherent state of lower polaritons. Note however that (cj^f-q)^ — 0, as c^, V-q 
are fermionic operators, thus this wavefunction can be also written as: 

I'^Bose) = CXp 

(Acos(fo)^i)n(l + ^ 

sin(Co)<^(q)cljU_q) |0) . (27) 

q 

Thus, if <f{q) has a step-like form, this can also describe a BCS-like state |96| . 
More generally, the parameters A, and the function (^(q) can be taken as variational 
parameters, and used to minimise the energy. 

Starting instead from the Hamiltonian of Eq. (n^) one is instead led to write: 



l^TLs) = e^^o Y[ (cos(^c«) + sm{9c.)Dl) {0)^^^ (28) 

a. 

= II^"(k)'^(q)cLek/(m,+m.)+q"m.k/(me+m.)-q (^9) 
k,q 

where we have now introduced $a(k) as the localised centre of mass wavefunction. 
Note that the operator, D^,, describing a localised exciton does not square to zero. 
It is thus not possible to rewrite the BCS-like product in Eq. (^) as an exponential; 
there is a qualitative difference between the states in Eq. (^7|) and Eq. ( ^8| ) . Although 
{■^aY 7^ 0> t;he product in Eq. (^) only allows each operator D]^ to occur at most 



Collective coherence in planar semiconductor microcavities 



18 



once, so for a given single-particle state labelled by a, only zero or one excitons may 
occupy it, and thus prevents multiple occupation. By including the disorder-localised 
centre of mass wavefunctions, Eq. ( |2^ ) describes single occupation of a set of localised 
exciton wavefunctions, while in comparison, Eq. (|2^) describes only the single, lowest 
energy, delocalised exciton mode. As above, we may take the parameters A, Oa and 
the functions $„(]<), (^(q) as variational. 

Unfortunately, direct evaluation of the expectation of Eq. (^) with these trial 
wavefunctions is challenging. At low enough densities, no multiple occupation occurs, 
so in this limit Eq. (^) and Eq. (^8|) become comparable: Expanding Eq. (^8|) for 
small 9a, the terms in the product can be rewritten approximately as: 



o{el). 



(30) 



This would be equivalent to Eq. ( Pq) except that Eq. (g6|) macroscopically occupies the 
k = exciton state, whereas Eq. (^0|) occupies a collection of disorder-localised states. 
Although not identical, a superposition of many localised states distributed across the 
sample can (at low enough densities) behave similarly to the translationally invariant 
k = state. Thus depending on the relative importance of disorder localisation, 
and on the difference of Coulomb interaction between different single-particle exciton 
states vs interaction for multiple occupation of the same single-particle state, one may 
find which of Eq. (^6|) or Eq. ( ^8| ) has lower energy. 

Furthermore, both of the above wavefunctions are mean-field approximations 
of the ground state, and in both cases, energy could be lowered by constructing 
the Nozieres-Bogoliubov state. To discuss this, let us consider the simpler case of 
structureless bosons, One can then understand this state in two ways, either as a 
variational ansatz, as in Ref. 



|A) = exp XL 



|0), 



(31) 



and then find A, Ak by minimisation. Alternatively, the same state can be described 
if one considers fluctuation corrections to the mean-field theory. As is well known, in 
the presence of a condensate, the quasi-particles are the Bogoliubov modes |Q, i.e. 

bI = cosh(<?!)k)i[ + sinh(<?!)k)i_k, (32) 
with (/)k the Bogoliubov rotation angle. Thus, given the Bogoliubov spectrum, 
the lowest energy state is the Bogoliubov vacuum, jfiBog), defined such that it is 
annihilated by all B^^, i.e. 

cosh(0k)-Lk + sinh(0k)-^L 



ll^Bog) = Vk 



which is clearly solved by: 



exp -Etanh(0k)4i4 |0> 



(33) 



(34) 



Two comments are in order about the significance of this state; firstly, the physical 
reason this state is of lower energy is the quartic interaction, in particular terms like 
L^L^_]^LqLq + H.C., which favour states which are not eigenstates of the number of 
k = particles. Secondly, even when projected to an overall number state, one may 
retain features of this state, by writing a superposition of terms with different division 
of the number of particles between the condensate mode and other states. 
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2. 3. Theoretical treatments — effects of the environment 

Having discussed various models of the polariton system, we now turn to how these 
models, and the effects of the environment, may be treated. We first briefly outline the 
thermal equilibrium case, and compare mean-field theories of the two models discussed 
above. We then discuss some of the various approaches that one may use to describe 
the effects of the environment, focusing mainly on non-thermal steady states. Finally, 
we try to separate and clarify the concepts of coherence, condensation, superfluidity 
and lasing, which while often related, need not necessarily occur together. 



2.3.1. Thermal equilibrium The simplest approximation for the environment is to 
consider the system in thermal and chemical equilibrium with a bath. While it is 
clear that the current experiments involve substantial pumping and decay, which will 
be discussed next, there are compelling reasons to deal with the equilibrium case. 
Firstly, the properties of a given model in the equilibrium case are instructive when 
considering the range of behaviour it can show; while the equilibrium properties of 
weakly interacting dilute Bose gas are well studied ||9^, the properties of 

models like Eq. (|l7|), with distributions of oscillator strengths and energies |7^, ^ are 
less known. Even within the weakly interacting Bose gas picture, interesting features 
can arise from considering non-quadratic dispersion | 10C1| , 101, 102|, or the effects of 
anisotropic spin interactions jij, |l03[ . The second reason is that with improvements in 
the quality of mirrors, and refinement to the design of microcavities and the conditions 
of pumping, experiments have been able to increase the thermalisation rate to be 
comparable to or faster than polariton decay rates and so for these, or future, 

experiments, the correct description may become increasingly close to equilibrium. 

The treatment of both Eq. (^) and Eq. ( |l7| ) in equilibrium can be put in a similar 
form by considering their saddle point, or minimum action equations. Formally, 
these can be derived by writing the imaginary time path integral for the partition 
function |104|, and then considering the configurations that minimise the imaginary 
time action. Thus, for the bosonic case within the effective lower polariton model, 
Eq. (|r^), the saddle point solutions satisfy a Gross-Pitaevskii equation: 



-idt + E, 



LP 
k=0 



pLP 



+ o(v^) + K:S,ol^ol 



(35) 



Here, as the dispersion El^ is not quadratic, we have expanded it to quadratic 
order to find the coefficient of V^. Note that by considering solutions of the form 



one can introduce the chemical potential, and thus recover the 



expected static Gross-Pitaevskii equation. 

For the fermionic model, more care is required; since there has been no projection 
onto lower polaritons, the saddle point condition leads to coupled_equation for the 
photon field and two-level systems 



-idt + ujQ ■ 



2m 



Using the spin notation of Eq. (jl^) one has: 
^o = E^^"'5(r-r„); 



dt^a — X Sq, 



Bq — 



.9q,o(V''o - 

gQ,oj(V'o 



-4) 



(36) 



(37) 
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Here is the localisation site of the two-level system S^. In the case where the only 
variation is ■00 (i) = ■0o6~*^* and the polarisation has the same time variation, one can 
eliminate the time variation by a gauge transformation. 

The sum over exciton energy levels can also be simplified if one makes two 
assumptions: firstly that the excitations are occupied according to a thermal 
distribution, and secondly that we can average over many realisations of excitonic 
disorder. This second assumption, that ■00 (r) varies slowly compared to the distance 
between excitons, or equivalently that the photon couples to many localised exciton 
modes, allows one to replace the sum over exciton energy levels with a sum over the 
statistical distribution of energies and the excitonic coupling strengths. This then 
yields: 



LOq- fl 



, V- .9a,otanh(/3£;„) 



2m ' ^ A 2Ea 



00 - , (38) 



where the energy = (ea — f^)'^ +ga olV^o(r)P depends on the local value of the slowly 
varying V'o(r). Note that in this way the exciton disorder does not lead to spatial 
inhomogeneity of the polariton condensate, and so in the absence of photonic disorder 
one would expect polariton condensation in the k = mode. This has a clear similarity 
to Eq. (|3^), but in this case, the non- linear interaction term is more complicated than 
it was in the bosonic case, V^°o ol-^oPj and the polariton-polariton interaction is due to 
the nonlinearity of the susceptibility arising from the saturable nature of the excitons. 
For a uniform and static condensed solution (V^'0o = = dtipo), the Gross-Pitaevskii 
equation (|3^) is also analogous to the gap equation (self-consistency condition) of the 
BCS theory pO^ . 

Despite their similarity, there is an important distinction between Eq. ( p5| ) and 
Eq. (38); Eq. ( |35| ) is temperature independent, while the nonlinear susceptibility in 
Eq. (3S) decreases at high temperature, and is eventually incapable of supporting 
condensation. Thus one can crudely say that Eq. ( |35| ) can support mean-field 
condensation at any temperature, and fluctuations must be considered to find 
a transition temperature. For the bosonic model, going beyond mean-field theory, 
one can also produce a temperature dependent equation by using a Hartree-Fock 
wavefunction, and thus including the effect of interactions between the condensate 
and non-condensed particles (see e.g. Rcf. for further details). In distinction 
Eq. (38) contains a finite mean-field transition temperature, and so fluctuations are 
only important when they significantly decrease this transition temperature |77|| . 
Thus, including fluctuations one finds a crossover from a fluctuation dominated phase 
boundary at low densities, to a phase boundary that is well described by mean- 
field theory in the high density limit, where long-range interactions dominate. The 
temperature dependence that appears in Eq. ( |38| ) can also be understood by noting 
that it was necessary to integrate out the excitonic degrees of freedom in Eq. (|3^) to 
produce an effective action for a single photon field. Thus, the saddle point density 
for the Bose-Fermi model contains both the condensate, and a contribution from 
incoherent excitons. 

By considering fluctuation corrections to the saddle point density, one may 
in three dimensions, and at low densities, recover the non-interacting transition 
temperature of a weakly interacting Bose gas. i.e. considering the transition from 
the normal side, a mean-field condensate density appears when /i — > 0, and so the 
critical temperature is given by ptotai = Pfiuct(7c,/i = 0). However, for a two- 
dimensional system, when the transition is of the Berezhinskii-Kosterlitz-Thouless 
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(BKT) class [105, 106, 107 , it is necessary to consider fluctuations in the presence 
of a quasi-condensate, as the BKT transition, where free vortices proliferate, occurs 
below the mean-field transition temperature. 

In order to calculate the total density in the presence of a condensate, it 
is important to note that the fluctuation corrections can deplete the condensate 
population, as well as increase the population of other modes [Q. The condensate 
density that comes from a mean-field calculation [i.e. from the uniform static solutions 
to Eq. (|35|)] is p = [io[^ = ^/V°^. One compact way of finding how fluctuations 
deplete the condensate density is by using the Hugenholtz-Pines relation as discussed 
in e.g. Ref. |9^ . Let us briefly summarise here how this argument shows that the 
condensate density is smaller than the mean-field estimate. To discuss this one must 
introduce the self energy of the condensate E. If we define the matrix of Green's 
functions: 

and introduce 5o(a;, k) as the free Green's function (i.e. in the absence of interactions), 
then the matrix of self energies is defined by E(w,/e) = QQ^{iu,'k) — Q~^{u!,'k). The 
Hugenholtz-Pines relation is the condition required of this self energy in order that 
there might be a gapless mode, as one expects for a Bose-condensed system. The 
condition can be written as Sn — E12 = fi. By writing the self energies in terms of 
the densities of condensate and normal state particles, one may use this identity to 
write the condensate density in terms of fi and the normal state density. At leading 
order in interaction strength, it can be shown [ p8| that the self energy is given by: 

En = 2F<=ff(po + Pi), S12 = V'^ipa - pi). (40) 
In this expression po is the condensate density, pi is the density of particles in 
other states, and pi is the anomalous density, X]k(-^k^-k)- ^^"^ thus finds: po = 
{n/V^^) — (2pi + pi), i.e. fluctuations reduce the condensate density below its mean- 
field value. 

To then extract the BKT temperature, one needs to find the condition for free 



vortices to proliferate [108, 107, 109|. As described in those works, this requires one to 
know the fugacity of a vortex, and the effective vortex- vortex interaction strength, both 
of which depend on the superfiuid stiffness p^, which may be found from the difference 
between transverse and longitudinal current response functions (see Eq. (|4|)). The 
result is a transition which occurs at fesTeKT = (2/7r) (p^/m). In the case of bosons 
with quadratic dispersion, and in the limit of weak quartic interaction strength V^^ , 
i.e. mV''^ < 1, one may extract an asymptotic relation between the superfiuid density 
and the total density, giving /cbTbkt = [2Tr/\n{B/mV°^)]{p/m), where quantum 



Monte Carlo calculations [11C| give B = 380 ± 3. The phase boundary calculated 
according to the fermionic model, i.e. using Eq . (|38| ), and the boundary for the BKT 



transition in a bosonic model following Refs. [IOC, 101, 102[, but with the effective 
incl usion of disorder, are shown in Fig. 0. These boundaries are discussed further in 



Sec. 2.4.1 



2.3.2. Pumping, decay, and non- equilibrium treatments A more realistic discussion 
of the experimental environment must consider that polaritons may escape, and 
so continuous pumping is required to produce a steady state. In addition, if one 
is to describe pulsed experiments, or the transient behaviour after the pump is 
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Figure 7. Phase boundaries for an equilibrium polariton condensate. Solid 
and dashed linea^ ruark the mean-field phase boundary of the boson-fermion 
model of Refs. [[73, M] (dashed line indicates the region in which this boundary 
is strongly modified by fluctuation corrections). Blue lines indicate the BKT 
transition temperature for thermal population of the lower polariton branch, 
with the modified Landau criterion as discussed in Ref. |102 . The interaction 



strength used for this curve is that which arises from the bosonic approximation 
to the boson-fermion model — i.e. includes the effect of disorder in the saturation 
interaction. For comparison, orange dashed line marks the phase boundary for the 
naive estimate (i.e. neglecting condensate depletion due to density fluctuations) 
BKT transition temperature for a quadratic dispersion with effective polariton 
mass. Plotted for Q.^ = 26meV exciton-photon detuning 5 = u)o — e* = 5.4meV, 



photon mass m 



ph 



2.59 X lO-^m- 



switched on, a dynamical approach is required to describe the time dependence of 

Considering for the moment steady state situations 



103|. 



population |111, 112, 

— i.e. c.w. (continuous wave) pumping — one may highlight two important features 
of the difference between the pumped, decaying system and thermal equilibrium. The 
first is that the distribution function; i.e. the population of each energy level, may be 
far from t hermal, and set inst ead by the balance of pumping, decay, and thermalisation 
rates; [El 



114 



115 , lis , 11£|. The second class of effects is that incoherent pumping 



and decay introduce dephasing, and can change the excitation spectrum of the system, 
the additional inclusion of these effects are discussed in Refs. |117, 116|, (see also the 
discussion in Sec. 3.2.l| ). There are a wide variety of approaches that may be applied 
to study one or either of these features; in the following we discuss briefly how some 
of these various approaches are related, and what limitations they may have. For a 



more general discussion see e.g. Refs. [118, pa, 119 



In order to describe the properties of the pumped, decaying system, one requires 
a method to calculate various correlation functions. Given an expression for single 
particle correlation functions, one may then find many properties of interest, for 
example the occupation of modes, the luminescence and absorption spectra, and the 
first-order coherence properties. The most general information about one particle 
correlations can be written in terms of the two correlation functions (t, r) = 
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('(/'U^i i")V'(0, 0)), ~ r)'0^(O, 0)), which with ij] describing the photon field 

correspond directly to luminescence and absorption probabilities. These encode 
information both about the form of the spectrum and thus the density of states, 
and also about the population of those states. For example, the density of states is 



given by Im[5 



[UJ ■ 



iO, k)], where the retarded Green's function can be written as: 



r) = Q{t)\Q^ (t, r) — Q"^ (t, r)]. In equilibrium, these two Greens functions can be 
related in terms of the thermal distribution function, but out of equilibrium no such 
simplification is possible. 

Let us now discuss different methods to calculate these Green's functions. The 
first method is to find and solve the operator equations of motion for "0, ^"^1 and thus 
to evaluate the correlation functions directly. In order to describe pumping and decay, 
one considers coupling the system to baths, which either pump particles and energy 
into the system, or provide modes into which particles may decay. These baths are 
assumed to be large, so their properties (e.g. distribution functions) are fixed, and 
not affected by the system. Since the bath and system are coupled, the equations of 
motion for the system operators will also include bath operators. If one considers the 
initial state of the bath to be drawn from some fixed (e.g. thermal) distribution, then 
the expectation of bath operators will be random quantities, with statistical properties 
set by the bath's distribution pSl; t hus such coupling to baths introduces noise, 
giving quantum Langevin equations |35|]. The second method is to write equations 
of motion for Green's functions, which will now involve correlation functions of bath 
operators, which again can be found if one assumes a fixed distribution for the baths. 
These coupled equations are the content of the Keldysh formalism |120|, in which 
it turns out to be simpler to write equations for various linear combinations of the 
Green's functions, allowing one to combine the Green's functions in a compact matrix 
notation (see e.g. |118|). The equations for the Green's functions can also be derived 
in a path integral formulation |119|. The path integral formalism allows one to make a 
close connection with the methods discussed in Sec. [2.3.1 , and can describe changes to 
both the occupation and to the spectrum induced by pumping and decay 1 117, 116|. As 
such, it allows one to introduce a complex self-consistency condition (i.e. a complex 
equivalent of the Gross-Pitaevskii equation); this interpolates between the laser, in 
which self-consistency requires balancing of pumping and decay, and the equilibrium 
condensate, where self-consistency instead relates real self-energy shifts. This point 



will be discussed in more detail in Sec. 3.2.1 



The above two approaches allow one to find self-consistently the population 
and the spectrum in the presence of pumping and decay. In certain cases — 
for example in the normal state, with weak pumping and decay — the changes 
to the spectrum may be small, and one is interested primarily in changes to the 
population. By assuming a known form for the retarded Green's function (i.e. for the 
spectrum of the system) it is possible to extract equations for the population from 
the Green's function formalism (see Ref. |118, Chapter 9] for details). These will lead 
to the Boltzmann equation [121|, which may also be derived phenomenologically, by 
considering various rates of transfer between different energies. Thus, by neglecting 
the "kinetic" effects, that change the spectrum, but retaining "dynamic" effects, that 
change the populations, one may investigate how pumping and decay change the 
occupation of the spectrum 114, 115, 112, 113|. Even without pumping and decay, 
the spectrum changes due to interactions in the condensed state [|l^; and further 
calculations retaining "kinetic" effects suggest the pumping and decay further modify 
the spectrum (see Refs. [117, 116| and also Ref. [122 in a different context). Because 



Collective coherence in planar semiconductor microcavities 



24 



changes to the spectrum modify the density of states, they can be expected to in turn 
affect the population dynamics; it is therefore not clear how valid it is to consider 
population dynamics in the condensed system without self-consistent treatment of the 
changes to the spectrum. 

A different approach to treating the coupling to baths is to consider the density 
matrix, which allows calculation of any single-time expectation of operators, as well 
as certain multi-time correlations, discussed further below. Calculating the full A^- 
body density matrix allows calculation of correlation functions of arbitrary numbers 
of particles, rather than just the single particle correlations discussed above. Density 
matrix methods, and Green's functions methods can be related as the single-particle 
density matrix (i.e. tracing over coordinates of all but one particle) is equivalent to 
the equal time part of the one particle Green's function. 

To numerically calculate evolution of the density matrix, due both to the system 
Hamiltonian, and to the effect of coupling to baths, one can choose an appropriate 
basis and write the density matrix in terms of a distribution function over this basis, 
i.e. p = ^ |a) (6| P{a,b). The time evolution of the density matrix then corresponds 



to the time evolution of this distribution function. Under certain conditions |35|, |123| , 
it is possible to interpret the evolution of this distribution as describing evolution of 
a quasi-probability distribution. In such a case, the equation of motion for P(a, b) 
is a Fokker-Planck equation, and can be rewritten in terms of Langevin equations 
for classical variables. If one chooses to resolve the density matrix onto a basis of 
coherent states, a variety of ways of doing this exist, among which we mention two 
important choices: the positive P distribution, and the Wigner distribution (See |Q 
for other possible choices, and further details). The positive P distribution formally has 
the desired properties, and gives the desired Langevin equations, but is numerically 
unstable when applied to the kind of problem discussed here [123|. The Wigner 
distribution, although not technically matching the above requirements, can have its 
equation of motion approximated — the truncated Wigner representation |123, 124| 
— which allows one to find appropriate Langevin equations, and is numerically stable. 

The Wigner distribution allows one to find the expectation of symmetrised 
products of operators; at equal times it is trivial to extend this to find general products 
of operators, by making use of equal time commutation relations. It is also possible 
to find multi-time correlation functions from such an approach, however since the 
unequal time commutation relations are not a-priori known, one cannot generally find 
expectations of other orders; in this sense the truncated Wigner approach does not 
allow G'^ and G'^, but only the symmetrised combination tj< +G'^ , and hence cannot 
separate the density of states from its occupation. 



2.3.3. basing, condensation, superfluidity Condensation, coherence and superfluidity 
are often connected, however when dealing with two-dimensional systems of composite 
particles, where finite-size and non-equilibrium effects may be important, it is 
important to separate and clarify these concepts. The discussion below compares 
condensation, coherence and superfluidity in equilibrium infinite systems to the case 
including finite-size and non-equilibrium effects. For simplicity we discuss these ideas 
in terms of a general bosonic field '0^, rather than any specific field appearing in the 
microcavity polariton system. The first concept is macroscopic occupation of single- 
particle wave- function; rigorously this can be defined as the existence of a macroscopic 
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eigenvalue of the reduced one-particle density matrix |125| 



pi(r,r') = (Vt(r)V'(r'))=E"^^*(r)¥'.(r'), (41) 

i 

where Ui is the occupation of the single-particle mode Lpi (r) . A macroscopic eigenvalue 
exists if \m\N^oonQ/N ^ where N = Y^^Ui is the total number of particles. In 
an infinite system, if the macroscopically occupied state is an extended state, then 
there is Off Diagonal Long Range Order — i.e. if in the position representation 
limr_,oo |iy9o(r)| 7^ 0, then there remain extensive terms far from the diagonal ]126| . 
Thus, in such a case, the correlation function limr^oo (V'U'")V'(0)) — "-o</'o('")'/'o(0) 7^ 
0. In a non-interacting two-dimensional trapped (and thus finite) gas of bosons 
there can be a sharp crossoverB to a state with a macroscopic eigenvalue (i.e. of 



the order of the number of particles) of the one-particle density matrix |127, 128|. 
However, the single-particle state (/3o(r) to which this eigenvalue corresponds will be 
a state localised in the trap, and so despite the existence of a macroscopic eigenvalue, 
limr^oo(V'^(r)V'(0)) =0. 

The visibility of interference fringes is directly related to the first-order coherence 
function gi{t = 0; r), where: 

, (■^'(-■')^(o.0)> 

^/(«.t(0,OW<(0,0))(*t(r,^W■(^,^)) 

coherence can be defined by the properties of this function. As just discussed, if 
one defines coherence by the limit of 51(0, r 00), then in a trapped system, this 
function vanishes. However, coherence will exist across the size of the trap. In such 
inhomogeneous and complicated binary classification of coherent / incoherent is 

less useful than a description of how the coherence varies as a function of separation in 
time and space. Both the size of this variation, and its functional form will depend on 
the interplay of finite size (and form of trapping potential), temperature, interactions. 



and decay rates |116| 



Superfluidity can meanwhile be defined separately as the diff erence of longitudinal 
and transverse response functions at vanishing wave- vector 1 12S , |l^ : 

X.,(c^ = 0,q) = 2 / dr((J,(q,T)J,(-q,0)» 
Jo 

= XT(9)(%-^)+XL(g)^, (43) 

and ps oc lim^_,o [XL{<f) — XT (9)]- Since superfluidity results from a change of the 
response functions, it occurs only in an interacting system; without interaction bosons 
do not become superfluid. In a two-dimensional infinite interacting system, below 



the BKT transition |105, 106, 107|, coherence decays as a power law rather than an 



exponential — low energy phase fluctuations prevent true long range order |130| 
and superfluidity exists, but no macroscopic occupation of a single mode. 

When one considers a more realistic system, which is both interacting, but also of 
finite extent, one cannot ignore a-priori the physics of the BKT transition, nor can one 

X In the limit of vanishing trap curvature and infinite number of particles the crossover becomes 



a phase trar 



particles [127 



tinr i. but only if trap curvature vanishes as the correct power of the number of 



128 



§ Above the BKT transition, in addition unbound vortices are present, and these lead to exponential 
decay of correlations. 
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ignore a-priori the physics of the trap. At low enough temperatures it is clear there will 
be macroscopic occupation of a single mode, and full coherence across the trap. How 
this state is approached as temperature is reduced, or as density is increased differs 
depending on whether interactions or finite size effects are dominant. If described as 
a non-interacting gas, the coherence at all distances increases uniformly as a single 
mode is increasingly occupied |127, 128|. In the BKT scenario, power law correlations 
develop on intermediate scales (between some short range thermal length and the trap 
size); then as temperature decreases, the thermal length increases and the p ower with 
which correlations decay decreases, again restoring full coherence as T ^ 131 132 1, 
as shown in Fig. @. 



Interacting, no trap 



Trap, no interactions 




0.5 1 
Distance 



Figure 8. Sketches of the spatial decay of coherence, g\{t = 0;r). Left; 
interacting infinite 2D Bose gas, showing power-law decay at long distances. 
Right: non-interacting trapped gas; at short distances (less than the thermal 
length) perfect coherence exists; on intermediate lengthscales coherence reaches 
a constant value, due to the non-zero condensate fraction; coherence eventually 
goes to zero for separations larger than the trap size. 



Adding non-equilibrium effects, the nature of decay of correlations in an infinite 
system is significantly altered |117, 116|. The long wavelength phase modes, 
responsible for decay of correlations become diffusive, i.e. the poles of the Green's 
functions, which in the equilibrium case have the form lo ~ ±cfc, take instead the 
form a; ~ ± {ck)'^ 



in the pumped and decaying case. This is illustrated 



in Fig. ||. Combining the effects of phase diffusion, and discrete level spacing [116| 
the properties of coherence are further modified, and one approaches the laser limit: 
temporal coherence of laser emission comes from one, or at most a few modes (resonant 
with the laser cavity), and so phase diffusi on of a s ingle mode leads always to 
exponential decay of temporal correlations |134, 135, 13£|. As well as this distinction 
of forms of temporal coherence, it is worth mentioning here a few other important 
distinctions between lasing and the generalised concept of condensation discussed here, 
as there are evidently also similarities |137]. Most obviously, polariton condensation is 
seen at the polariton resonance, which is significantly (of the order of the Rabi splitting 
Q,r) below the lowest cavity photon mode; as such nonlinear emission coexisting 
with strong coupling is a signal that one should consider polaritons, and not just 
photon lasing. A more fundamental difference is that lasing requires inversion, while 
condensation does not; this is a consequence of the standard laser systems possessing 
little coherence in the gain medium, while excitons, being part of a coherent polaritons, 
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Figure 9. Sketch of momentum dependence of real and imaginary parts of 
poles of the Green's function, in the presence of pumping and decay. At zero 
momentum, there is a free phase, i.e. a mode with vanishing real and imaginary 
part. For small momentum both modes are diffusive, with no real part, but a 



non-zero imaginary part. Only above some critical momentum d oes 
develop. Such results are seen both in the call 
also for resonant pumping (as discussed in Sec. 



ons in R n fs 1 
2.l|) in Refs. [ 12^ ' 



L33 . 



lie] 



,1 part 
and 



are coherent |138, 139, 140, 141, 142]. This distinction can also be seen by comparing 
the critical lasing condition to the Gross-Pitaevskii equations of Eq. ( |35| ) and Eq. (p8|). 
In the presence of pumping and decay, the susceptibility (describing the nonlinear 
response of excitons) becomes complex [117, 116 ; the real part of the susceptibility 
gives the nonlinearity in the Gross-Pitaevskii equation. In contrast, the imaginary part 
of the susceptibility describes absorption or gain, and leads to the lasing condition, 
that round trip gain and loss balance. A treatment of a model system with pumping 
and decay elegantly shows how these conditions can be combined, giving an expression 



in terms of the total susceptibility [117, 116[. 

Starting from the strong-coupling regime, when crossing over to the weak-coupling 
regime (see, e.g. [143, gM 144|), the polariton splitting collapses, and so the lasing 



mode no longer has any excitonic character, and becomes the standard photon laser. 
A major advantage of an exciton-polariton laser over standard lasers is that it can 
operate without the inversion of the electronic population [145|, and therefore it has 
a much smaller threshold pump power. It is interesting to note that wide band gap 
semiconductors, such as GaN and ZnO, would be particularly suitable as in these 
cases excitons are stable at higher temperatures and densities, and therefore they 
could operate at room temperature. Electronic population inversion is not necessary 
because in the exciton-polariton laser, both photon field and excitons are coherent. 
In addition, the involvement of the excitonic field leads to strong nonlinear effects 
compared to conventional lasers, due to exciton-exciton interactions. 



2-4-. Phenomena 

The previous two sections discussed the models, and the treatments of the 
environment, that have been used to theoretically model polariton condensates. This 
section in contrast will review a few of the phenomena that have been predicted as 
possible signatures and properties of a polariton condensate. 
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2.4- 1- Tc and the phase boundary Within a given model, and effective description 
of the environment, it is natural to first ask under what conditions a condensate can 
exist. Within an equilibrium model of the lower polariton branch as weakly interacting 
bosons, phase diagrams for the physical parameters of various possible materials are 
shown in Refs. |10C, 101| (however, see Ref. |102| for a discussion of the effects of non- 
quadratic dispersion on the BKT transition temperature) . By considering a simplified 
version of the Bose- Fermi model Eq. (p^), where the energies are described by a 
Gaussian distribution, while all excitons display a fixed coupling to light, the mean- 
field phase boundary was first calculated in Refs. js^, The effect of fluctuations, 
restoring the bosonic limit at low densities was instead considered in Refs. jo^, |77|] . 
Since the content of the boson-fermion model at small densities and temperatures is 
equivalent to a bosonic model, and the low momentum part of the polariton dispersion 
is controlled by the photon mass, it is not surprising that it is possible to recover the 
standard BKT transition temperature of a weakly interacting Bose gas from the boson- 
fermion model in the low density limit. A calculation of the mean- field boundary which 
instead takes into account a realistic description of the quantum well disorder and the 
full distribution of oscillator strengths (see Fig. |^) has been performed in Refs. j?^, Q 
(see Fig. 0). 

Owing to finite size effects the experimental systems do not have a sharp phase 
transition marking the onset of a broken symmetry. All the observed transitions are 
rounded, and in order to extract a phase boundary from experiment, some criterion 
has to be chosen. One commonly used criterion is the nonlinear threshold; i.e. the 
point at which the relation between emission at k = and input pump power becomes 
nonlinear. Such a criterion is somewhat problematic. A second-order phase transition 
can be expected to be accompanied by a region with large susceptibilities, and thus 
such nonlinearity extends over a significant range of parameters, and so identification 
of a strict phase boundary from it is hard. Only in a mean-field theory does the 
onset of nonlinearity occur at the transition. However, because of the long-range 
nature of interactions, mean-field theory can be an adequate description for lasers, and 
for polariton condensates except at very low densities (as can be seen from Fig. 
The following sections discuss other phenomena that may demonstrate or describe 
condensation and coherence in microcavity polariton systems, and may thus provide 
alternative, or corroborating experimental criteria to find the phase boundary. 



2.4-2. Energy-resolved luminescence, resonant Rayleigh scattering In both the weakly 
interacting boson model, and the Bose-Fermi model, condensation leads to changes 
in the spectrum of polariton modes — most significantly the appearance of the phase 
modes In addition, for a model starting from localised exciton states, with 

a distribution of energies and oscillator strengths, there is weak emission from sub- 
radiant exciton states between the upper and lower polaritons |146|| . This emission is 
also modified by condensation [^: in the presence of a condensate these sub-radiant 
exciton states have energies Ea as defined following Eq. (^), and so the density of 
states is changed by the coupling to the coherent field. This change to the emission 
is discussed further in detail Ref. |Q. In practice, it is however hard to observe 
the incoherent luminescence of thermally excited modes in the presence of a strong 
signal from the coherent condensate. One suggestion to overcome this problem is to 
probe these excited modes via resonant Rayleigh scattering; by using a phase sensitive 
measurement one may be able to identify a small coherent scattering signal even in 
the presence of emission from the condensate Fig. |lO| shows the Rayleigh 
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Figure 10. Contourplot of the disorder averaged RRS intensity (/pq(i^)) for 
p| = |q| as a function of the dimensionless momentum |p|aox and rescaled 
energy 2{u) — /i)/f2j^, for zero detuning, Rabi spUtting Oj^ = 26meV, temperature 
ksT = 20K, and a disorder strength characterised by an inverse scattering 
time 1/t = 1.16meV. (a) non-condensed (density p ~ 0); (b) condensed regime 
(p ~ 3.6 X lO^cm"^) ; (c) condensed regime (p ~ 3.1 X lO^^cm"^) . Green lines 
mark chemical potential: if non-condensed RRS emission is present only above 
the chemical potential, but when condensed, it is present both above and below. 
[From [£§]. 



scattering signal expected both above and below the phase transition; note in the 
condensed case, one sees linear modes both above and below the chemical potential; 
this is as one expects from the Bogoliubov spectrum, where the normal modes are 
superpositions of particle creation and annihilation. 



2.4-3. Momentum distribution of radiation By integrating the luminescence at a 
given wave-vector, one may consider the momentum distribution of polaritons — 
this has also been discussed in the context of exciton condensation [132, 147 1. In a 
weakly interacting, two-dimensional, infinite system, the density of states would be 
constant, and this would just show the Bose-Einstein distribution. In experiments, the 
distribution at small momentum is close to, but can deviate from the Bose-Einstein 
distribution However, such deviations are expected; both due to interactions, 

which modify the spectrum and so modify the density of states ^ , and also due 
to finite size, which cuts off components at small wave-vectors. A naive calculation 
shows how the change of spectrum causes a change of density of states: for a weakly 



Collective coherence in planar semiconductor microcavities 



30 



interacting Bose gas with dispersion e\^, and mean-field equation /i — gpo, we have the 
Bogohubov spectrum = ± \/£k(£k + '^gpo), and simple algebra gives the density of 
states: 



z/s(tj,k) = Im[g^(cj + iO+,k)] =fm 
-d{uj ~ Ek) H 



iLu + ek + gpo 

uj^ + El 
Ek- ei^- gpo 



S{uj + E,,) .(44) 



2-Ek 2£'k 
This suggests that as k 0, the density of states goes like gpo/E]^ cx 1/fc, however this 
neglects the fact that the low energy modes are phase modes, and phase fluctuations 
may grow without bounds, as only their gradient costs energy; a full calculation, 
e.g. 1 141, 77 , shows that this term becomes l/fc^"'' with ry ~ 2T/Tbkt- In order to 
include the effects of finite size, one approach is to start with the zero-temperature 
Thomas-Fermi spatial profile, for which coherence across the whole cloud leads to 
sharp angular peaks, and then consider how phase fluctuations destroy the long range 
coherence, and thus soften the peaks |132|. There have also been calculations of this 
momentum distribution for non-equilibrium situations, where details of the dynamics 
of polariton relaxation lead to a maximum of the distribution at non-zero k, e.g. 
Refs. m fll |l|, pT2|. 



2.4-4- Linewidth, first- and second-order coherence, polarisation The line-shape of 
emission is controlled by the Fourier transform of the first-order coherence function 
as a function of time, i.e. gi{t\Y = 0), using giit^r) defined in Eq. (p2[). A perfectly 
coherent single mode source would have gi{t\Q) constant, and thus infinitely sharp 
lines. Although at any finite temperature, population of slow phase modes will lead 
to some decay of gi{t]Q), the transition to a condensed phase will lead to a slower 
decay of 3i(i;0)|j], as has been seen Similarly, as discussed in Sec. 2.3.2, one can 
also consider the spatial decay of coherence ; whic h ha s also been much studied in 
cold atomic gases (see, theoretical discussion in Refs. |148, 131, 149 1 and experiments 



in Refs. |150|, |l5lj, |l5| ) 



The calculation of coherence has already been discussed briefly in Sec. 2.3.3. As 



was stressed there, a distinction exists bet ween cohe rence of a few mode laser [134], 
applied to the polariton problem in Refs. |135, 136 1; and coherence in a continuum 
of modes 111, |ll7|, In both cases increasing temporal coherence in the system 



should lead to a narrowing of linewidth as the phase transition is approached, however 
the behaviour far above the transition may differ. This is clear from the fact that for 
a single mode with a quartic interaction, the equilibrium state is a number state, and 
so starting from a coherent state, one has phase diffusion due to self phase modulation 
even in the absence of any pumping or decay [ |136[ . For a many-mode system, the 
ground state is neither a coherent nor a number state, and is better described by the 
Bogoliubov-Nozieres state |^^, thus it is not clear that the same effect — i.e. self 
phase modulation causing larger broadening — should persist. 

A related measurement is the second-order coherence function (see e.g. Ref. Q): 

(Vnr,0)^t(r,t)^(r,t)V;(r,0)) 



52 (t) 



(V't(r,0)^(r,0)) (V't(r,i)^(r,t)) 



(45) 



II One should note in defining a coherence time that, for an equilibrium, infinite, two-dimensional 
system, the long time decay would be power law, and so there is no well defined coherence time, 
however finite size effects lead to exponential decay 116|. 
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As mentioned above, for a thermal state, 52 = 0) = 2, while for a coherent state 
52(i = 0) = 1. Experimental measurement of 52 — 0) is restricted by finite detector 
integration times, and since 32 (i) = 1 at times when gi{t\Q) —f 0, it is hard to 
distinguish the value of 52 = 0) |3^, The dynamical behaviour of 52 — 0) 
as a function of time following switching on the pump laser ]153| , 154, |155| , or from 
an equilibrium picture with separate coherent and incoherent contributions has 
been considered. Within the Bose-Fermi model, the second-order coherence has been 
studied for a finite number of excitons coupled to a single coherent field, wh ere the 
finite number of states replaces the phase transition with a smooth crossover [ 156 |. 

Considering the differences of interaction strength between polaritons with 
parallel and anti-parallel spin polarisations, Laussy et al. |Q have shown that 
condensation should be associated with spontaneous development of a linear 
polarisation. That a condensate should choose some definite polarisation state is 
a much more general result, i.e. that repulsive interactions prevent fragmentation of a 
condensate, even when there are two single-particle states with identical energies |]96|| . 
Considering the specific form of the interaction, Laussy et al. showed that this 
specific polarisation should be a linear state (as opposed to circular or elliptical). 
Note, however, that the pinning of the polarisation to one of the crystallographic axes 
observed in Refs. j4|, ^ ^ cannot be an effect of spontaneous symmetry breaking, 
but it is instead likely due to some optical anisotropy of the microcavities. In atomic 
gases, investigation of "spinor condensates" is hard, as it requires the use of all optical 
trapping, as a magnetic field would Zeeman split the atomic spin states; further the 
populations of spin species are effectively fixed at the start of the experiment, so rather 
than polarisation, phase separation into spin domains has been observed there [157 . 



3. Resonant pumping 



Since its first realisation in 2000 |158|, the possibility of reaching the stimulated 
scattering regime for polaritons by resonant pumping has attracted considerable 
interest [|l59[ |l6^, |l6^, |l6|, |l| |l6|, |l6|, |l6|, [l66[ ^ |l^, |l6|, , and has initiated 
the search for polariton lasers and a new generation of ultrafast optical amplifiers and 
switches. In resonant excitation experiments, polaritons are optically pumped at an 
energy and momentum which allows coherent polariton-polariton scattering directly 
into the ground state. In pulsed, ultrafast pump-probe experiments, a second pump 
is used to initiate the stimulated scattering process to the ground state, while for 
c.w. (continuous wave) excitation the stimulated regime can be self-induced above 
some pump threshold. As in non-resonant excitation experiments, the direct mapping 
between in-plane momentum and angle of bulk photons is crucial. It is this direct 
mapping that allows one to directly pump at a given momentum, and thus to perform 
the experiments described here. 

The macroscopic occupation of the ground state in these experiments differs 
in some ways from condensation and spontaneous coherence arising from incoherent 
pumping. In experiments without a probe pulse, there is a free phase between the 
pump and signal modes, however as discussed in Sec. 3.2.1 , the existence of a free 
phase alone is not sufficient to ensure superfiuidity in such driven systems. Further, 
since the coherence of the signal beam is directly coupled to the pump beam, higher 
order correlations, and the linewidth of the signal may be inherited from the pump. 
Experiments with an additional probe beam at the signal frequency differ yet further; 
in such experiments the phase of the signal is fixed by the probe, and so in this case 
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Figure 11. 

^2k, 



(a) Lower polariton dispersion; (b) counterplot of l-E^^ + 
k — 2_Bkjj I as a function of k and the zero value contour (white 
dashed). A pair of final states (signal and idler) can be found by intersecting 
the white dashed curv e with a straight line passing through kp (indicated 
as pump). [From [17S] Copyright (2001) by the American Physical Society]. 



all excitations ought to be gapped. 

It is interesting to note that the idea of parametric scattering discussed here for 
polaritons has been recently a ppli e d in a different field; that of dilute gases of atoms 
confined in an optical lattice |l71 



later in Sec. 



4.1 



172 1 . We will discuss this example in more detail 



3.1. Summary of experiments 

The main idea of resonant pumping experiments |15S, 159, |16(]| , 161| , 162, 162, 164, 
|165| , |166| , |167] , |16§| , |16g| , |17C| ] is that of the coherent scattering of two polaritons from 
the resonantly pumped mode (pump) into the ground state (signal) and a high energy 
state (idler) [see Fig. |ll|(a)]. Energy and momentum conservation in this scattering 
requires one to have {kp,kp} i— )■ {0,2kp}, where 

E]^l,+E\l^2E\f^ , (46) 

which uniquely selects the momentum (i.e. angle) of pump, signal and idler. If one 
instead relaxes the condition k = for the signal mode, then for a fixed pump angle kp, 



it can be shown 173 , 168 1 that those final states which satisfy energy and momentum 
conservation describe a figure-of-eight in momentum space [see Fig. pl](b)]. The 
non-parabolic dispersion of polaritons is crucial in order to obtain such a resonant 
scattering processes, and so the analogous process for excitons is forbidden. Resonant 
experiments can be divided into two types: In the first type the scattering is stimulated 
by a weak probe field (parametric amplification) , while in the second type there is no 
probe and the seed to initiate stimulation is provided by the pump itself, if above a 
certain power threshold (parametric oscillation). 
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Figure 12. Contour maps of the PL emission spectra for excitation at the 
'magic' angle (x) (a) when no probe is applied at zero angle and (b) when a 
probe is applied, which stimulates a strong gain in the probe direction. The 
peak of the PL (white line) shows a marked asymmetry and a blue-shift 
of ~ 0. 7me V with respect to the low power polariton dispersion (dashed). 
[From 1 159 1 Copyright (2000) by the American Physical Society]. 



In their pioneering work [158|, Savvidis and collaborators realised the stimulated 
scattering regime for the first time in an InGaAs/GaAs/AlGaAs microcavity with a 
Rabi splitting oiQ,B. = 7meV. By pumping polaritons at the "magic" angle ~ 17° (for 
zero exciton-photon detuning) close to the inflection point of the LP dispersion curve, 
and using a second pump at zero angle (probe) to initiate the process, a substantial 
gain of up to 70 was observed. The large observed signal required both the "magic" 
angle pumping, and the probe at zero angle, and was absent either if no probe was 
applied ]159| (see Fig. or when pumping at different angles even in presence of 
the zero angle probe. These results provide strong evidence for a polariton scattering 
process stimulated by the probe. The bosonic scattering rate is enhanced by a factor 
N ^ 1 where N is the population of bosons in the final state |174|, in this case the 
ground state. Moreover, stimulated scattering was confirmed in this experiment by 
the observed exponential dependence of the gain on pump power. Note however that 
the obtained stimulated regime is more correctly described as a 'parametric' amplifier 
rather than a laser, i.e. the population of the lowest mode is amplified by a phase 
coherent parametric process |163|. 



161, 162 



Much experimental work has followed this first result |159, 160 
163, 164, 165, 16(:, 167, 168, 169, 1701. Evidence of stimulated scattering has been 
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obtained in plj |, with a three beam pulsed experiment, where two pumps excite states 
at large and opposite angles, 9^ ^ ±45°, and a third beam is used as a probe at 
normal incidence. Here, the pump creates a quasi-thermal exciton reservoir at large 
momentum; then, scattering of polariton from kp and — kp to the k = upper and 
lower polariton states can occur, and is is stimulated by the occupation (due to the 
probe) of the final state. Note that, at zero detuning, such a scattering process 
conserves energy {2e'^ — Eq^ + Eq^) and momentum. Recent experiments in the 
original two beam pump probe configuration, but with time resolved measurements 
as a function of different pump and probe angles have beg un the in vestigation of the 

hydrodynamic properties of the injected "polariton fluid" [17£, 17£]. 

In another series of experiments, making use of res onant c.w. excitation 162 160 | 



rather than ultrafast pulsed excitation as in |158, 159, |3l[, the stimulated scattering 



regime has been reached even without the probe beam (parametric oscillation). Here, 
the stimulated scattering is self-initiated when the pump power is strong enough that 
the final state population becomes close to one. It is interesting to note that, as in 
the non-resonant excitation experiments, a minimum pump power is required (with 
threshold-like behaviour) in order to overcome the 'bottleneck' effect and thus provide 
the occupation of the ground state required to stimulate scattering. In these particular 
experiments, occupation of the ground state was estimated to be close to 1 at threshold 
and around 300 at the highest pump powers. A maximum blue-shift of the lower 
polariton of the order of 0.5meV was observed, which is much less than the Rabi 
splitting of ~ 6meV, and so confirms that the experiment is always in the strong 
coupling regime. Because of the stimulated scattering, pumping of the ground state 
mode is efficient, and so these systems have a low power threshold, typically 5 times 
smaller than for a high quality vertical cavity surfaces emitting laser (VCSEL) |160|. 

While the experiments described above were at temperatures of the order of 
5K, Saba and collaborators have reported pump- probe parametric amplification of 
polaritons with an extraordinary gain up to 5000 at temperatures up to 120K in 



GaAlAs-based microcavities, and up to 220K in CdTe-based microcavities [164 . The 



highest possible operating temperatures for observing amplification have been shown 
to be determined by the exciton binding energy (25meV for the used CdTe wells and 
13.5meV for GaAs wells), rather than the polariton Rabi splitting, which in these 
experiments varies from 25meV for the 24 quantum-well (QW) CdTe microcavity to 
15meV and 20meV for GaAlAs microcavities with respectively 12 QWs and 36QWs. 
Time-resolved measurements, obtained by controlling the delay between pump and 
probe, show ultrafast dynamics of the parametric gain, promising future applications 
in high-repetition-rate optical switches and amplifiers. 

As an alternative to the finite angle resonant configuration, a few experi- 
ments [163, 169 have concentrated on the double energy and momentum resonance. 



where pump, signal and idler are all at normal incidence, kp — 0. In this geometry 
it is possible to investigate the dependence of the phase of the signal on the phase 
of the pump, and thus to show [ [163| that in this degenerate configuration the para- 
metric scattering is a coherent process. In addition, in most of the early resonant 



experiments, such as Rcfs. [159, |l60|, 162, 164|, spectral narrowing of the signal was 



observed when amphfication occurs, also suggesting that the signal is coherent [166|. 
Direct evidence for the coherent nature of the signal emission has come only recently, 
in Ref. |170[ , where first-order temporal and spatial coherence were investigated. In 
that work, in a resonant c.w. pump, no probe, configuration, two spots separated by 
70/zm, coming from the same laser excited region of the sample (lOO^m FWHM of 
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Gaussian profile), are overlapped in momentum space, showing interference fringes. 
In addition, by making use of noise measurements, the emitted signal is shown to be 
in a single-mode quantum state, rather than in a multi-mode state. 

Recently, the pair correlation of the emitted signal-idler polaritons have been 
demonstrated by showing that polaritons in two distinct idler modes can interfere if 



and only if they share the same signal mode |177 



3. 2. Theories 

Predictions about threshold conditions, spectral properties and efficiency of the 
amplification in resonant pump-probe experiments can be easily obtained by making 



use of the effective lower polariton Hamiltonian i/LP ( MD described in Sec. 2.2.2, to 
which one must add the coupling to the external radiation pump r2pump(0 and probe 



^^probc(i) fields |17§: 



i^cxt = I [^k,kpf^pmnp(t) + f^k.of^probcW] COS 6'ki[ "f H.C.j . (47) 

k 

A closed set of equations of motion for the expectation values of probe (signal) 
{Lo{t)), pump ^Lkj,(i)), and idler (L2kp(0) niodes can be obtained by factorising 
field expectation values, and neglecting higher order correlations. By solving such 
equations, both in the s tead y-state regime and numerically for the pulsed excitation. 



Ciuti and collaborators [178| have obtained the conditions for gain threshold, showing 
that the efficiency of the amplifier depends very strongly on the polariton linewidth; 
the larger the linewidth, the higher the threshold and the lower the maximum gain. 
Similar results, such as the blue-shift of the signal with increasing pump power. 



have been obtained by Whittaker [179| in a classical nonlinear optics treatment. In 
Whittaker's paper, a phenomenological model — where a nonlinear excitonic oscillator 
is coupled to the cavity mode, driven by external fields — can describe both parametric 
amplification and oscillation. This treatment, considering classical fields, and retaining 
only the frequencies corresponding to pump, signal and idler modes, is equivalent to 
factorising the field expectation values and considering the equation of motion for 
(Lo(t)), (Lk,(i)), and (L2k,(i)). 

By considering the case of a c.w. pump without a probe, and expanding up to 



the second-order in the field expectation values, it can be also shown [173, uM that 



below threshold for parametric amplification, the polariton photoluminescence, 

/•oo 

PL(k,t,w) ~ cos2 6'k3?y rfTe-^("-*^)^^L[(< + T)Lk(t)) , (48) 



has a blue-shifted asymmetric emission distribution, as was observed in |15£], and 
shown in Fig. The lower polariton blue-shift here. Elf E^f + Vk,kp,o (^kp)^ 
is due to the interaction term of Eq. (|l6|). Note that the above formalism is valid 
only below the threshold for parametric emission, as above threshold the equation of 
motion for the pump mode, describing pump depletion should also be included. 



Recent theoretical work [180 , including the effects of polariton blue-shifts on the 
parametric oscillator equations, has shown that in the c.w. configuration the 'magic' 
angle is not necessary, as it is in the ultrafast pump-probe case, and that, under 
suitable pump conditions, the parametric oscillator can in general be observed for 
pump angles 9 > 10°. Such calculations seem to explain the experimental results 



obtained in Ref. [167[. In that experiment, with c.w. pumping, stimulation was 
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achieved over a wide range of pump angles, from 10° to 24° and with the signal always 
at k = 0. Including the polariton blue-shift also allows one to distinguish parametric 
oscillation, where output intensity grows continuously, and bistability, where output 
suddenly jumps at some threshold pump power, and shows hysteresis if power is then 
reduced. The regimes of bistability and parametric oscillation have been investigated 
as a function of pump angle for resonant pumping |180[ , and also allowing for mismatch 



between pumping frequency and the polariton energy at the pumping angle |181| 



There has also been a proposal to study the hydrodynamic properties of the 



injected polariton fluid by studying its coherent scattering by disorder | 182 |; i.e. 
resonant Rayleigh scattering. In the proposed experiment, a strong pump beam creates 
a large coherent population of polaritons, and also provides the source of polaritons 
which may be coherently scattered: In the presence of disorder, the polaritons can 
resonantly scatter to states with different momentum, but the same energy. This can 
be observed by looking at the angular distribution of photons escaping the cavity 
which are resonant with the pump beam. At low pump power, this emission pattern 
will be a ring, at the pump angle — i.e. those states with the same energy have 
the same modulus of momentum. However, since a large coherent population of 
polaritons modifies the polariton dispersion relation, at higher pump powers, both 
the shape and intensity variation across the pattern of resonantly scattered photons 
reveal information about the polaritons in the cavity. 

3.2.1. Phase degree of freedom and low energy modes The laser and an equilibrium 
polariton condensate form extreme ends of the spectrum of systems in which 
coherent emission results from a symmetry breaking transition; the resonantly pumped 



polariton laser falls somewhere in between. Although, as discussed in Sec. 3.1, the 
coherence of the signal is inherited from the pump, in the parametric oscillator 
configuration, without a seed signal beam there is in principle a free phase A0 between 
the signal and pump modes. That is, the equations of motion are invariant under the 
transformation Lq L^e^^'^, L2i<i — > L2W e~'^^'^, using the notation introduced on 
page In an equilibrium condensate, the invariance of the energy under a global 
rotation of phase implies the existence of a soft phase mode. However, the existence 



of a free phase does not necessarily lead to the same consequences 122 117 , 133 116 1, 
as we will discuss next. 

Amongst several distinctions between a laser and an equilibrium condensate, one 
important difference is that for a laser, the threshold condition is the balance of gain 



and decay 137, 174 1, while for a condensate, the Gross-Pitaevskii equation (gap 
equation, self-consistency condition) involves the real parts of the self-energy. This 
difference implies a difference in the spectrum of fluctuations. Both cases have a pole in 
the fluctuation spectrum at w = and zero wavevector, corresponding to global phase 
rotations. However, for small, but non-zero wavevector, the spectrum first acquires a 
real part for a true condensate, describing linearly dispersing modes |128| , p8| (in an 
interacting system); if instead one must balance gain and decay, the spectrum instead 



first acquires an imaginary part, describing diffusive modes [122, 117, 133, 116 (see 
Fig. Although both the laser and strongly pumped polariton systems share this 
diffusive structure, there is an important difference in the properties of the phase 
mode between a conventional laser and the strongly-pumped polariton system. The 
difference is that, as discussed in Sec. |2.3.3 , a conventional laser has a discrete spectrum 



of wavevectors (all other modes rapidly decay), while the microcavity polariton system 
has a continuum of in-plane wavevectors with comparable decay rates. The continuum 
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of diffusive phase modes should thus lead to differences between the coherence 



properties of a laser and a resonantly pumped polariton system ]122| , 133]. The 
difference between diffusive and dispersive low energy degrees of freedom may also 
have implications for pattern formation in nonlinear systems [ 185| , 184]. 



3.2.2. Polarisation and spin relaxation By considering the spin of polaritons, with 
resonant pumping, one can consider the coupled dynamics of the polarisations of the 
pump, signal and idler modes. This dynamics leads to a rich variety of physical effects, 
due to the interplay between spin dependent stimulated scattering, and precession 
induced both by momentum dependent TE-TM (Transverse-Electric, Transverse- 
Magnetic) splitting, and other energy splittings due to polariton-polariton interactions. 
The following is a brief summary of the theories that have been applied to explain 



these fe atur es. Interest in the subject began with experiments, in both c.w. Jl85] and 



pulsed ]186] experiments, that c ould not be explained by regarding each spin species 



as acting independently. In Ref. ] 185 ], signal intensity as a function of pump ellipticity 
(from linear to circular) was studied, and a maximum found at an intermediate 
value. In Ref. |l86| , a large output signal was seen for a linearly polarised pump, 
but a circularly polarised probe, and was explained there in terms of stimulated spin 
scattering. In addition, the direction of the linear component of polarisation of the 
signal was observed to vary as a function of the degree of ellipticity of the pump 



beam. This large rotation was discussed in Ref. ]187] as a giant Faraday effect, with 
the cavity amplifying the effect of spin splitting of the cxciton energy levels, with 
the spin splitting being due to unequal spin populations of the exciton states. In 



contrast, Ref. ]188] showed that parts of the above results could be explained by 
introducing coupling between the two circular polarisations. Such coupling provides 
two new terms; parametric scattering of cross polarisations, and processes where pump 
and signal, or pump and idler polaritons exchange polarisation. These terms lead 
to a threshold power that depends on ellipticity of the pump, and can under some 
conditions show lowest threshold for an elliptic polarisation. 

In a different experiment, Martin et al. ] [189| investigated the dependence of the 
circular polarisation of the non-linear emission on the detuning between the cavity 
photon modes and excitons. This result, showing oscillations of the polarisation, have 



been attributed to the dependence of the TE-TM splitting on the detuning ] 19C ] 



This TE-TM splitting, which depends on wave-vector, can lead to time-dependent 



oscillations of the polarisation of the signal and idler modes, as discussed in 1191 
In order to combine the precession due to such a splitting with a spin-dependent 
stimulated scattering to the signal state, one is led to write a spin dependent 



Boltzmann equation ]191]. This formalism was sufficient to reproduce the results 



of the experiment described in Ref. ]192]. In that experiment the circular polarisation 
component of the signal was studied as a function of pumping strength, for both linear 
and circular polarisation. It was found that, for a linear pump, there is a maximum 
of the degree of circular signal polarisation near the nonlinear threshold. As discussed 



in Ref. ]192], this effect can be understood as a competition between two effects; self- 
induced Larmor precession, which rotates the pseudospin describing the polarisation 
from linear to circular polarisation, and stimulated scattering to the ground state: Far 
above threshold, the rate of scattering to the ground state is too high to allow time 
for any change of polarisation in the pump mode. 

An ingredient missing from the works listed so far was stimulated scattering 
due to polariton-polariton interactions (as opposed to polariton-phonon scattering. 



Collective coherence in planar semiconductor microcavities 



38 



stimulated by final state polariton population). Shelykh et al. considered the 



dominant, parallel spin, interaction in Ref. ]193{. However it was later realised [194 



195 1 that the scattering of anti-parallel spin states, though small, is important, as 
it leads (at high densities) to a 90° rotation of linear polarisation direction between 
pump and signal. In a recent experiment the polarisation of the output signal 
was studied for a linearly polarised pump, as a function of pump power and angle of 
the linear polarisation. This experiment showed that, just above threshold, the signal 
was elliptical, with the major axis (the direction of the linear component) rotated 
with respect to the linear polarisation of the pump — i.e. the direction of major axis 
depended on the direction of polarisation of the pump, but the angle between the two 
directions was not constant. In addition, the degree of linear polarisation of the signal 
decrease d far above threshold (as well as the degree of circular polarisation as observed 



in Ref. 192 |). This rotation was reproduced in a model combining precession due 
to static and self-induced (Larmor) splittings (including an extra in-plane splitting, 
as discussed in Ref. |]67| ) as well as spin rotation in stimulated polariton-polariton 
scattering. The reduction of linear polarisation far above threshold was explained by 
rapid self-induced Larmor precession, which rotates the linear polarisation direction 
so rapidly that it averages to zero. 

4. Connection to other systems, conclusions 

^.1. Atomic gases and Feshbach resonances 

Microcavity polaritons and their condensation are related to the physics of two- 
component atomic Fermi gases near Feshbach resonances. In particular, the crossover 
from a fluctuation dominated phase boundary to a mean-field phase boundary with 
increasing polariton density is closely related to the BEC-BCS crossover recently 



studied in these atomic gases [ 196 , 197 [. For atomic gases, in contrast to polariton 
experiments, the density of particles is typically kept fixed, while the interaction 
strength is varied via magnetically tunable Feshbach resonances, allowing one to go 
from a BEC of bound molecules to a BCS state of fermionic pairs. The interaction 
strength is tuned by changing the detuning between the zero of energy for pairs of 
atoms in their original spin states, and a closed channel resonance level of atoms in 



some higher energy spin configuration (See Ref [198[ and Refs therein for more details). 
The Bose-Fermi model used for polaritons and described in section 2.2.3| is very 



similar to the model initially proposed for the description of the BEC-BCS crossover 



in atomic Fermi gases [199[. The fermionic operators h\a'^ in the polariton model, 
introduced in Eq. (p^), are analogous to the two spin species of atoms in the Feshbach 
resonance model; the photon is analogous to the closed channel resonance level; 
the dipole interaction relates to the hyperfine interaction; and the polariton to the 
Feshbach molecule. There are, however, a number of important differences between 
the two models; most notably the absence of a direct four-fermion interaction in the 
polariton model, and the existence of an energy dependence, and a distribution, of 
exciton-photon coupling strengths. Secondly there is a marked difference of mass 
ratios: For polaritons, the photon mass is typically a factor of 10^ times smaller 
than the exciton mass; in the Feshbach case the closed channel resonance has a 
dispersion controlled by a mass which is twice the atomic mass. Note also that there 
is a difference in interpretation between the photon and the closed channel resonance 
level in Feshbach resonance. Although inside the microcavity a photon is not an 



Collective coherence in planar semiconductor microcavities 



39 



eigenstate — the upper and lower polaritons are instead eigenstates — outside the 
cavity it is possible to physically separate the photons from excitons. Similarly, the 
closed channel resonance level is not an eigenstate; however, in distinction from the 
photon, it cannot be physically separated from the other two-body states to which it 
is coupled: In general one cannot have a given molecular resonance level in isolation 
from other molecular states. Another difference between the systems comes in how the 
gap equation [i.e. Eq. (|38|)], and constraint on the density can be used to find chemical 
potential and temperature: In the mean-field (BCS) limit of the atomic case, the closed 
channel resonance level lies at high energies, the chemical potential lies well within the 
band of fcrmionic states, and so the density depends only on chemical potential, and 
not on temperature. In the polariton model, the chemical potential remains below 
the bottom of the band of fermionic states, and so both temperature and chemical 
potential influence the density. Rather, the fluctuation dominated (BEC) to mean- 
field (BCS) crossover of the polariton model is in the nature of which excitations are 
responsible for depopulating the condensate at finite temperatures. 

Another issue worth considering is that further analysis has questioned the need 
to use the Bose-Fermi model in applications to experimentally relevant Feshbach 
resonances. It has been shown |200| , ^Ol| that the molecules created in the vicinity 
of a Feshbach resonance are halo dimers extending over large distances in which 
the closed channel admixture is tiny. Thus, the resonance level acts to enhance 
the effective interaction between fermionic pairs, but the crossover to loosely bound 
molecules does not rely on the macroscopic occupation of this level (as was initially 
suggested). Since the resonance level is not significantly occupied, it was shown. 



by using realistic atomic potentials [202|, that the BCS-BEC crossover in atomic 



Fermi gases near Feshbach resonances is of the same nature as originally considered 



by Leggett |203]. That is, it is based on the smooth crossover of the pair size and not 
on the macroscopic occupation of the resonance level. Thus, microcavity polariton 
experiments in the normal operating conditions of large photon field occupation 
present the first experimental realisation of the BEC-BCS crossover which differs 
substantially from the original scenario p03| p|[77|. 



Finally, it is interesting to note that the same idea of parametric scattering and 
amplification in a resonant pumping configuration as is described in Sec. ^, has been 
recently applied to an ultracold dilute gas of bosonic atoms confined in an optical 



lattice ||I72| (see also the related e xper iments on dynamic instabilities [204, 205 1, and 
four wave mixing of matter waves |l7l| ). There, a ^^Rb Bose- Einstein condensate was 
loaded into a moving one-dimensional optical lattice. The optical lattice causes the 
atomic dispersion to deviate from quadratic, and allows parametric scattering: Atom 
pairs with initial momentum kp inherited from the moving lattice scatter elastically 
into two final states; kp — k and kp + k. By generating a seed of atoms with momentum 
kp — k, parametric amplification of both the seed and the conjugate momentum has 
been observed, with a gain determined by the atomic scattering length. 

4-2. Excitons, quantum hall bilayers, triplons 

The high quantum degeneracy temperatures and the high degree of control obtained 
by laser photo-excitation suggest that excitonic systems should also provide excellent 
environments in which to study macroscopic coherence phenomena. Much interest 
over the last two decades has been attracted by excitons in bulk CU2O crystals (see. 



e.g.,Ref. \20t\ and references therein). Excitons in CU2O have a large binding energy 
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(~ 150meV) and the fact that the direct dipole transition from the exciton ground 
state is forbidden guarantees low radiative recombination rates. In CU2O crystals 
however non-radiative recombination processes such as the Auger effect cause loss and 
heating, and represent the main obstacle to the observation of quantum degeneracy 
in these structures. 

Low radiative recombination rates together with high cooling rates can be 
obtained for spatially indirect excitons in coupled quantum wells. Here, an electric 
field is applied along the growth direction, in such a way that electrons and holes 
separate into different wells. In contrast to the case of direct excitons, lifetimes up to 
few microseconds can be achieved, while the high cooling rate gives a much shorter 
thermalisation time, typically in the nanosecond range. As a consequence, in coupled 
quantum well structures, thermal equilibrium with the lattice can be relatively easily 
obtained by either waiting a few nanoseconds after photo-excitation or allowing the 
excitons to travel away from the excitation spot. While much experimental work has 
been done on these structures in the last few years [207, 208 , 209 210| (for review see 
[^ll| , |212| ), unambiguous evidence for Bose-Einstein condensation of indirect excitons 
in coupled quantum wells is still missing and is the subject of intensive on-going 
studies. Recently, gases of indirect excitons have been trapped and equilibrated 
using in- plan e potentials, either by applying localised str ess t o change the local band 
energies [ [213[ or by means of optical |214] or electrostatic |215| traps. The confinement 
of indirect excitons prevents the fast reduction of initial density which occurs in the 
absence of trapping due to fast expansion driven by their strong dipole-dipole repulsive 
interaction and their relatively high mobility. 

It is interesting to note that, for an untrapped system, while the formation of an 
external excitonic ring has been explained as an in-plane charge separation, at low 
temperatures it has been shown that, such a ring can separate into a periodic array of 
beads, and that the light emitted by each of these beads is coherent ]21(| |. The origin 
of such a phenomenon is still unknown. 

Recently, experimental evidence suggests that, under appropriate conditions, 
an electron-electron semiconductor bilayer system in the Quantum Hall regime can 
condense into a superfluid state which might be interpreted as an excitonic-like 
superfluid |217|. In a bilayer 2D electron system with total filling factor vt — 1, 
excitonic pairs can be thought of as formed by filled electron states in one layer and 
empty electron states in the second layer. By changing the ratio between intra- and 
inter-layer Coulomb interactions, signatures of the transition to a condensed excitonic 
phase have been shown by a dramatic increase in the tunnelling rate between the two 
layers at zero interlayer voltage |218|, in Coulomb drag measurements [219 , and in 



counterflow measurements 220 



Another condensed matter system in which a phase transition can be described as 
condensation of an excitonic mode is that of magnetic "triplon" excitations. This has 
been seen in a variety of compounds including TlCuCla [221, 222 , BaCuSi206 [223[ 



and Pb2V30 [224[, where by changing the applied magnetic field, there is a crossing 
between spin singlet and spin triplet excitations. The resultant magnetic phase 
transition can be described as condensation of the triplet magnon mode. Very recently 
such a transition was seen at room temperature, using parametric laser pu mpin g to 
create a non-equihbrium density of triplons in the compound Y3Fe2(Fe04)3 |225 . 
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4-3. Conclusions 

In this article we reviewed the experimental results to date demonstrating coherence in 
microcavities, and discussed the variety of theoretical models and techniques that have 
been used to describe it. We discussed experiments with both non-resonant pumping, 
in which coherence may spontaneously arise from an initially incoherent source of 
polaritons, and the optical parametric amplifier and optical parametric oscillator: 
Both kinds of experiments allow one to explore the interplay of strong-coupling, 
coherence, lasing and condensation. In our discussion of theoretical descriptions, 
we highlighted those aspects of these solid state systems which can introduce new 
questions about coherence — disorder, decoherence, particle flux, potentially non- 
thermal distributions. 

Of course, there are subjects connected to coherence in microcavities that we have 
not had space to discuss, or have only discussed briefly. For example, we have only 
briefly mentioned here questions about the dynamics of condensate formation, and of 
the polariton response following short pump pulses. The last years have seen rapid 
experimental progress in this field, with the first convincing evidence of coherence 
developing from incoherently injected polaritons in a variety of systems pl| , ^ |58| . 
This experimental progress both gives hope for the possibility of future experiments, 
and applications on coherent microcavity polaritons, as well as focusing attention on 
those areas in which further theoretical work is necessary. There are areas in which our 
discussion has been brief because some questions have only been partially addressed 
to date: Questions about hydrodynamics in such partially coherent, pumped decaying 
systems. In order to address questions about the generic behaviour of such systems, 
it is important to understand how the variety of models used to describe polaritons 
relate, and what features each can explain. Finally, let us mention that there remains 
an interesting topic which can be the subject of much future research: what new 
experiments are possible in these light-matter systems that were not possible in either 
lasers or in atomic gases. 
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